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ABSTRACT 

We present extensive ugrizY H JKg photometry and optical spectroscopy of 
SN 2005gj obtained by the SDSS-II and CSP Supernova Projects, which give 
excellent coverage during the first 150 days after the time of explosion. These 
data show that SN 2005gj is the second clear case, after SN 2002ic, of a ther- 
monuclear explosion in a dense circumstellar environment. Both the presence of 
singly and doubly ionized iron-peak elements (Fe 111 and weak S 11, Si 11) near 
maximum light as well as the spectral evolution show that SN 2002ic-like events 
are Type la explosions. Independent evidence comes from the exponential decay 
in luminosity of SN 2005gj, pointing to an exponential density distribution of 
the ejecta. The interaction of the supernova ejecta with the dense circumstellar 
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medium is stronger than in SN 2002ic: (1) the supernova hnes are weaker; (2) 
the Balmer emission hnes are more luminous; and (3) the bolometric luminosity 
is higher close to maximum light. The velocity evolution of the Ha components 
suggest that the CSM around SN 2005gj is clumpy and it has a flatter density dis- 
tribution compared with the steady wind solution, in agreement with SN 2002ic. 
An early X-ray observation with Chandra gives an upper-limit on the mass loss 
rate from the companion of M < 2 x 10^* Mq yr~^. 

Subject headings: supernovae: general — supernovae: individual (SN 2005gj) 



1. Introduction 



Thermonuclear supernova explosions (Type la supernovae, SN la hereafter) are believed 
to be the detonat ion or deflagration of a white dwarf accreting matter from a companion 
star (lArnettl Il982l ). The mass of the white dwarf slowly increases until it approaches the 
Chandrasekhar limit where the star becomes thermally unstable. At this point fusion of 
Carbon and Oxygen begins near the center and quickly moves through most of the star before 
degeneracy is lifted. The result is a spectacular and powerful explosion that is visible across 
much of the Universe. Since SN la arise from a narrow range of white dwarf masses, th eir peak 



lumin osities are very consistent and they make excellent distance indicators (e.g., iPhillips 
19931 ). SN la are powerful probes of cosmology and have been instrumental in narrowing the 
uncertainty in the Hu bble parameter, disc o yery of the accelerating u n iverse , and constraining 



dark energy models (IHamuv et al. 



1999; Riess et al. 200 



i 



Astier et al. 



1995 



1996 



i 



Riess et al 



2006; Riess et al 



2006 



2005 



1998 



Per. 



Wood-Vasey et al 



mutter et al. 



20071). 



But the use of SN la as reliable distance indicators will always be questioned until the 
progenitor and explosion physics are well-understood. What types of binaries create SN la? 
How is matter transferred to the white dwarf without causing thermonuclear runaways on 
the surface? Are there several types of progenitors?. These big questions remain to be 
answered and detailed observations of hundreds of events have yielded few clues. 



In 2002, iHamuy et al.l (120031 ) identifled a new kind of supernov a. The early spe ctrum 
of SN 2002ic was a cross between a Type la event and a Type Iln fiDeng et al.ll2004j . have 
called this type of supernova a "Ha"), showing P-Cygni features similar to SN la and resolved 
Balmer lines in emission. Ty pe Iln supernovae are core-collapse explosi ons going off in dense 
circumstellar environments (jSchlegellll990l : IChevalier fc Franssonlll994j ). They are relatively 
common since the massive stars that create core-collapse supernovae often have thick winds. 
If the interpretations of the pre-explosion observations of SN 2005gl are correct, SN Iln could 
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be associated in some cases with luminous blue variables (iGal-Yam et al.ll2006l ). 



In the case of SN 2002ic, the presence of Balmer lines with profiles characteristic of 
SN Iln, and the high luminosities and slow decline after maximum lead to the conclusion 
that most of the energy came from the interaction of th e ejecta with a dense circumstel- 
lar medium (GSM ). Other Type Iln events (SN 1997cy, ICermanv et"aD I2OO0I : SN 1999E, 



Rigon et a. 



2003r) have been re-classified as SN 2002i c-like that were caught late in their 



evolution (IHamuv et al.ll2003l : IWood-Vasev et all 120041 ) . 



SN 2002ic provided the first direct evidence that thermonuclear explosions can also oc- 
cur in a dense medium, but in t his case the circumstellar mediura is pr o bably generated by an 
Asymptotic Giant companion (IHamuy et al.l l2003l : IWang et al.l 12004 : iHan &: Podsiadlowski 
20061). Howe v er, th ere is still debate in the literature about the origin of SN 2002ic. 
Livio &: Riesd (120031 ) proposed the merger of two white d warfs as a po s sible progenitor 



with the explosion occurring in the common envelope phase. iGhugai et al.l (120041 ) concluded 
that the properties of the circur nstellar interaction in 2002ic-like events can be broadly ex- 
plained by the SN 1.5 scenario (jiben &: Renzinil 119831): the thermonu clear explosion of the 
degenerate core of a massive AGB star. Recently, iBenetti et al.l (120061 ) questioned the earlier 
interpretation of the observations and proposed that SN 2002ic can be equally well explained 
by the core collapse of a stripped-envelope massive star in a dense medium. 

SN 200 5gi was discovered on 2005 September 28.6 (UT) by the SDSS-II Supernova Gol- 
laboration (iFrieman et al.l 120071 ) in gri images obtain ed with the SDSS-2.5 m telescope at 
Apache Point Observatory (APO). The new supernova (IBarentine et al.ll2005l ) was ~1" from 
the center of its host galaxy at the position a = 03^^01"^12!0, S = +00°33'13'.'9 (J2000.0). It 
had SDSS magnitudes {g,r,i) = (18.6, 18.6, 18.7) mag, obtained from PSF photometry after 
kernel-matching and subtraction of a template image in each band. The S N was indepen- 



dentl y discovered by the Nearby Supernova Factory on 2005 September 29 (lAldering et al. 
2006h . 



SN 2005gj was classified as a Type la candidate from the first three epochs of the gri 
light curve using a light curve fitting program, and was sent to the queue of the MDM-2.4m 
telescope for spectroscopic confirmation. The optical spectrum obtained on 2005 October 
1 (UT) showed a blue continuum with resolved Hydrogen Balmer lines in emission, very 
similar to the spectrum of a young Type Iln supernova, but also with an unusual continuum 
showing broad and weak absorption features. Further spectroscopic follow-up showed a 
dramatic evolution. The continuum became substantially redder and developed broad, P- 
Gygni features probably associated with blended lines of Fe-peak mass elements, similar to a 
Type la SN a few weeks after maximum. The spectrum obtained on 2005 Nov. 12 (UT) was 
remarkably similar to that of the unusual Type la supernova SN 2002ic obtained on 2002 
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Dec. 27 (UT) flPrieto et al.lboosh . 



Aldering et al.l (120061 ) presented optical photometry and spectroscopy of this SN. Through 
detailed analysis they confirmed its photometric and spectroscopic resemblance to SN 2002ic, 
confirming it was a new case of a Type la explosion interacting with a dense circumstellar 
environment. From a spectrum obtained with the slit oriented to overlap with its host galaxy, 
they calculated a redshift for the host o f z = 0.0616 ± 0.0002. S N 2005gj was not detected 
i n the radio with th e Very Large Array (jSoderberg &: Frail! 120051 ) or in the X-ray with Swift 
Jimmler et allboosl ). 



Here we present extensive follow-up photometry and spectroscopy of the Type la/Type Iln 
SN 2005gj during the first ~150 days after discovery. These are the most detailed observa- 
tions ever obtained of a SN 2002ic-like event, and provide insight into the early evolution, 
progenitor, and variety of these events. We also present a sensitive, early X-ray observation 
with Chandra that gives an upper limit on the X-ray luminosity of this peculiar object. We 
describe the optical and NIR photometry of SN 2005gj in § [2] and the optical spectroscopy 
in § [31 We describe the X-ray observation with Chandra in § HI An analysis of the photo- 
metric and spectroscopic data are presented in § O Finally, we discuss the significance of 
the results in § [61 We will present later photometric follow-up and Spitzer/IRAC observa- 
tions in Paper II. We adopt a cosm ology with Hn = 72 ± 8 km s~^ Mpc~\ = 0.3, and 
Qa = 0.7 throughout the paper (see lSpergel et al.ll2006l ). which yields a distance modulus of 
fi = 37.15 mag to the host of SN 2005gj. 



2. Photometry 



2.1. SDSS and MDM 



The Sloan Digita.1 Sky Survey (SDSS; [York et al. 



200oh uses a wid e-field, 2.5-meter 



telescope (iGunn et al.l l2006l ) and mosaic CCD camera (iGunn et al.l Il998l ) at APO to sur- 
vey the sky. The SDSS-II Supernova Survey, part of a 3-year extension of the original 
SDSS, uses the APO-2.5m telescope to detect and measure light-curves for a large num- 
ber of supernovae through repeat scans of the SDSS Southern equatorial stripe (about 
2.5 deg wide by ~120 deg long) over the course of three 3-month campaigns (Sept-Nov. 
2005-2007). SN 2005gj was discovered in the second month of the first campaign (Octo- 
ber 2005). Twenty epochs of ugriz photometry were obtained between 2005 Sep 26- Nov 30 
(U. T. ) . Details of the photom et ric system, magnitu d e system, astromet r y and calibra t ion ar e 
given inlpu 



Pier et al. 



cugita et al. 



Jl996h . iLupton et alj I Jl999h.[Hog"g et al.l J200ll ). ISmith et all J2002h 



J2003h . Ilvezic et al.l J2004I ). and 



Tucker et al 



J2006h . 



Additional griz imaging of 
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SN 2005gj was obtained with the MDM Obse rvatory 2.4m telescope using a facihty CCD 
imager (RETROCAM; see lMorgan et al.ll2005l . for a complete description of the imager). 



Photometry of SN 2005gj on the SDSS ir nages was carried out u sing the scene modeling 
code developed for SDSS-II as described in iHoltzman et al. (120071) . A sequence of stars 



around the supernova were taken from the list of llvezic et al.l (120071 ). who derived standard 
SDSS magnitudes from multiple observations taken during the main SDSS survey under 
photometric conditions. Using these stars, frame scalings and astrometric solutions were 
derived for each of the supernova frames, as well as for the twenty five pre-supernova frames 
taken as part of either the main SDSS survey or the SN survey. Finally, the entire stack of 
frames was simultaneously fit for a single supernova position, a fixed galaxy background in 
each filter (characterized by a grid of galaxy intensities), and the supernova brightness in 
each frame. 

The supernova photometry in the MDM frames was also determined using the scene 
modeling code. Since the MDM observations had different response functions from the 
standard SDSS bandpasses, the photometric frame solutions included color terms from the 
SDSS standard magnitudes. To prevent uncertainties in the frame parameters and color 
terms from possibly corrupting the galaxy model (here affecting the SDSS photometry), 
the MDM data were not included in the galaxy determination, but the galaxy model as 
determined from the SDSS was used (with color terms) to subtract the galaxy in the MDM 
frames. The resulting SN photometry from the MDM frames is reported on the native MDM 
system, since the color terms derived from stars are likely not to apply to the spectrum of 
the supernova. 

Figure [1] shows a 3.5' x 3.5' field around SN 2005gj and 16 comparison stars used for 
calibration of the SN magnitudes by SDSS and the Carnegie Supernova Group (CSP; see 
§ 12.21 below for details). In Table [T] we present the SDSS ugriz and CSP u'g'r'i' photometry 
of the comparison stars in common. The final SDSS and MDM griz photometry is given in 
Table [2 



2.2. CSP 

Optical photometry from the CSP was obtained with the Swope-lm telescope at LCO, 
using a SITe CCD and a set of u'g'r'i' filters. A subraster of 1200 x 1200 pixels was read 
from the center of the CCD, which, at a scale of 0.435 "pixel"^, yielded a field of view 
of 8.7' X 8.7'. Typical image quality ranged between 1" and 2" (FWHM). A photometric 
sequence of comparison stars in the SN field was calibrated with the Swope telescope from 
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observations of SDSS standard stars from lSmith et al.l (120021 ) during four photometric nights. 
SN magnitudes in the u'g'r'i' system were obtained differentially relative to the comparison 
stars using PSF photometry. In order to minimize the contamination from the host galaxy 
light in the SN magnitudes, PSF-matched ugri template images from SDSS were subtracted 
from the SN images. On every galaxy-subtr acted image, a PSF was fitted to the SN and 
comparison stars within a radius of 3". See iHamuy et al.l (120061 ) for further details about 
the measurement procedures. 

The NIR photometry of SN 2005gj was obtained by the CSP using three different instru- 
ments/telescopes. A total of 15 epochs in F, J and H filters were obtained using RetroCam, 
mounted on the Swope-lm telescope at LCO. A few additional ep ochs in YJHKg were ob- 
tained with the Wide Infrared Camera ( WIRC: iPersson et al.ll2002[) mounted on the duPont- 
2.5m telescope, and the PANIC camera (IMartini et al.ll2004l) mou nted on the Magellan- 6. 5n i 



Baade telescope, both at LCO. We refer to iHamuy et al.l (120061 ) and iPhillips et al.l (120071 ) 
for details of the imagers and procedures to extract the SN photometry. The host galaxy 
was not subtracted from the NIR frames; therefore the SN photometry contains an unknown 
galaxy contamination component. 

The final CSP u'g'r'i' photometry of SN 2005gj is given in Table [3] and YJHKg pho- 
tometry in Table HI A minimum uncertainty of 0.015 mag in the optical bandpasses and 
0.02 mag in the NIR is assumed for a single measurement based on the typi cal scatter in 



the t ransformation from instrumental to standard magnitudes of bright stars (IHamuy et al. 
2006h . 



3. Spectroscopy 

The spectroscopic observations of SN 2005gj are summarized in Table [61 They were 
obtained using five different telescopes/instruments at four observatories. 

A total of twelve spectra were obtained between early-October 2005 and late- January 
2006 with the Boiler & Chivens CCD Spectrograph (CCDS) mounted at the MDM-2.4m tele- 
scope. This instrument uses a Loral 1200x800 pixel CCD with 15 /xmpixel"^ and a 150 1/mm 
grating (blazed at 4700 A). We used a 2" slit which gives a dispersion of 3.1 Apixel"^ in the 
wavelength range ~3800-7300 A. 

Eight spectra were obtained at LCO with the Wide Field Reimaging Camera (WFCCD) 
and the Modular Spectrograph (ModSpec) at the duPont-2.5m telescope, and LDSS-3 at the 
Clay-6.5m telescope. In the case of WFCCD a 400 1/mm grism and a 1.6" slit were used, 
reaching a dispersion of 3 Apixel"^ in the wavelength range ^3800-9200 A. For ModSpec we 
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used a 300 1/mm grating and a 1" slit that gave a dispersion of 2.45 A pixel ^ in the range 



3800-7300 A. Further details of WFCCD and ModSpec can be found in iHamuv et al. 



((20061). For LDSS-3, which employs a STA0500A 4064 x 4064 pixel CCD, we used the VPH 
blue and red grisms, the latter with an OG590 filter to block second order contamination, 
with a 0.75" slit reaching a dispersion of 0.70 Apixel"^ and 1.12 Apixel"^ in the blue and 
red sides of the spectrum, respectively. 

Additionally we obtained two early spectra with the Double Imaging Spectrograph (DIS) 
mounted on the ARC-3.5m telescope at APO, and one spectrum with the Intermediate 
dispersion Spectrograph and Imaging System (ISIS) at the WHT-4.2m telescope at the 
Roque de Los Muchachos Observatory in La Palma, Spain. The DIS spectrograph has blue 
and red detectors, each uses a Marconi 2048 x 1024 pixel CCD with 13.5 /impixel"^. We 
used a 300 1/mm grating and a 1.5" slit which gives a dispersion of 2.4 Apixel"^. The ISIS 
spectrograph has a blue and red arm, the blue arm using a EEV12 CCD and the red using 
a MARC0NI2 detector. We used a 300 1/mm grating in both the blue and red arm and a 
1" slit, which gives a dispersion of 0.86 Apixel"^ in the blue and 1.47 Apixel"^ in the red. 

Most of the spectra were obtained close to the parallactic angle to minimize relative 
changes in the calibration of the blue and red parts of the spectrum due to differential refrac- 
tion through the atmosphere. The spectroscopic reductions were performed using standard 
IRAF tasks and included: bias and overscan subtraction, flat-fielding, combination of 2-4 
individual 2D spectra to reach the best signal-to-noise ratio in the final image, tracing and 
extraction of a ID spectrum from the combined 2D image, subtracting the background sky 
around the selected aperture, wavelength calibration using an arc-lamp, and flux calibration. 
In order to flux calibrate the spectra we observed 1-2 spectrophotometric standard stars per 
night. The spectra from LCO, WHT and MDM were corrected by atmospheric telluric lines 
using the spectrum of a hot spectrophotometric standard star and the spectra from APO 
were corrected using a model atmosphere. This correction is not optimal for some of the 
spectra and there are evident residuals left in the corrected spectra. 

Figure [2] and [3] show a montage of the optical spectra of SN 2005gj obtained from 
October 2005 to March 2006. We have split them in two figures to avoid crowding. The 
position of the most conspicuous spectral features have been indicated in this figure. 



4. X-ray Observation 



SN 2005gj was observed under Director's Discretionary Time for 49.5 ks on 2005 Dec 
11/12 (ObsID 7241) with the Chandra X-ray Observatory's Advanced CCD Imaging Spec- 
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trometer (ACIS). The data were taken in timed-exposure mode with an integration time of 
3.2 s per frame, and the telescope aimpoint was on the back-side illuminated S3 chip. The 
data were telemetered to the ground in "very faint" mode. 

Data reduction was performed using the CIAO 3.3 software provided by the Chandra X- 
ray Centei0. The data were reprocessed using the CALDB 3.2.2 set of calibration files (gain 
maps, quantum efficiency, quantum efficiency uniformity, effective area) including a new bad 
pixel list made with the acis_run_hotpix tool. The reprocessing was done without including 
the pixel randomization that is added during standard processing. This omission slightly 
improves the point spread function. The data were filtered using the standard ASCA grades 
(0, 2, 3, 4, and 6) and excluding both bad pixels and software-flagged cosmic ray events. 
Intervals of strong background flaring were searched for, but none were found. 

Absolute Chandra astrometry is typically good to ~0'.'5, and we sought to improve it 
by registering the Chandra image with an SDSS image. Chandra point sources were found 
using the wavdetect tool, and their positions were reflned using ACIS Extract version 3.101. 
Fourteen X-ray sources had SDSS counterparts, which we used to shift the Chandra frame 
by a small amount (O'.'IS in RA and O'.'OT in DEC). After the shift, the residual differences 
between the Chandra and SDSS sources had rms values of 0'.'19 in RA and 0'.'12 in DEC. 

We extracted counts in the 0.5-8 keV bandpass from the position of the supernova using 
a standard extraction region (~1" radius), and we constructed response flies with the CAIO 
tools and ACIS Extract. The background region is a source-free annulus centered on the 
position of the supernova with inner and outer radii of 6" and 32". Based on the 300 photons 
detected in this region, we expect 0.3 background counts in our source extraction region. 

We detect only two counts from the location of the supernova, but neither may be 
associated with the supernova itself. The counts had energies of 4.0 keV and 6.5 keV, but 
one would expect some emission in the 0.5-2.5 keV range since this is where Chandra has 
the most collecting area. For example, the average effective area in each of the 0.5-2.5 keV, 
4.0-6.0 keV, and 6.0-8.0 keV bands is 470 cm^, 270 cm^, and 90 cm^, respectively. 



We calculate our upper limits using the Bayesian method of iKraft et al.l (Il99ll ). For the 
0.5-4 keV band, the 68% (95.5%) upper limit to the source counts is 1.14 (3.05). For the 
0.5-8 keV band, the 68% (95.5%) upper limit is 3.52 (6.14). 

Since no Type la supernova has been conclusively detected in X-ray, we have no a priori 
expectation of the spectral shape. We therefore adopt a simple absorbed power law with a 
photon index of 2 and an absorbing column of uh = 7.08 x 10^" cm~^. For this choice of 



^http : //asc . harvard . edu 
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spectrum, the count rate to flux conversion is 5.2 x 10 erg/count in the 0.5-4 keV band 
and 6.9 x 10^^^ erg/count in the 0.5-8 keV band. 

We therefore arrive at 68% (95.5%) upper hmits on the X-ray luminosity of 1.1 (2.9) x 
]^q39 gj^g g-i ^Yie 0.5-4 keV band and 4.3 (7.6) x 10^^ erg s~^ in the 0.5-8 keV band. Based 
on the above statements concerning the Chandra effective area, we feel the 0.5-4 keV limit 
is a more appropriate limit. 



5. Results 

5.1. Optical light curves and colors 

Figure m shows the early SDSS and MDM ugriz light curves combined with the late time 
coverage given by CSP u' g'r'i'Y JHKs photometry. They give excellent multi-wavelength 
optical and NIR coverage and sampling of the first ~150 days after discovery (~140 rest- 
frame days after the time of explosion). 

In Table [5] we give important parameters derived from the light curves. We have a good 
estimate of the time of explosion of the SN at JD 2,453,637.93±2.02 (September 24.4 UT, 
2005) calculated as the average time between the epoch when the supernova is detected at 
> 5(7 in all filters at JD 2,453,639.94, and the last pre-discovery observation of the field 
at JD 2,453,635.91. The times and observed magnitudes at maximum in different filters 
are presented in Table O They are calculated from a high order polynomial fit to each 
light curve. To estimate the errors we assume a gaussian distribution for the magnitude 
uncertainty at each epoch and filter (assuming they are not correlated from epoch to epoch), 
then we draw randomly ~1000 simulated light curves and fit each one with the same high 
order polynomial. The la uncertainties are taken as the rms deviation of the mean values 
calculated from each simulated light curve. 

The risetimes, defined as the time between explosion and maximum light, become longer 
at redder wavelengths. They are 13.5, 19.7, 33.7, and 46.9 days in ugri filters, respectively, 
which correspond to 12.7, 18.5, 31.7, and 44.2 days in the rest-frame of the supernova. 

We also give in Table [5] the rest-frame magnitudes at maximum in different filters. They 
have b een corrected by Ga lactic exti nction in the l i ne of sight, using E{B — V)Gai = 0.121 



(SFD: ISchleeel et all Il998r ) and the ICardeUi et al.l (119891 ) (CCM) reddening law assuming 



Ry = 3.1, and i^-corrections (see below), which are not negligible in this obj ect since it is at 



redshift z ~ 0.062. We assume that the host galaxy extinction is negligible (lAldering et al. 



20061). 
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i^'-corrections have been calculated from the multi-epoch spectra presented in § [31 In 
order to estimate accurate /^-corrections we need good spectrophotometric calibration in all 
the wavelength range. Figure shows the differences between the observed g — r photometric 
colors obtained from the light curves and the synthetic colors calculated directly from the 
spectra, as a function of the observed photometric colors. We use an 8th order polynomial 
fit of the light curves to obtain the observed g — r color at the epoch of a given spectrum to 
better than ~ 0.03 mag. We can see that most of the spectra have good spectrophotometric 
calibration in the wavelength range oi g ,r filters (3800-7000 A) , with a residual color of 
~0.05 mag rms (difference between observed and synthetic colors), but there are some ob- 
vious outliers. To correct the spectra and produce a better spectrophotometric calibration 
consistent with the observed spectral energy distribution (SED) obtained from broadband 
photometry we warp the spectra multiplying by a smooth function to match the observed col- 



ors, a technique commonly used for calculating i^-corrections in Type la SN (INugent et al. 



2002|). 



First, we extrapolate the continuum in the blue and red sides of each spectrum presented 
in § |3] with a low order polynomial to have complete coverage of the SDSS ugriz filters 
(2000-11000 A). Using this extended version of the spectra we apply the CCM reddening 
law iteratively until the synthetic g — r color matches the observed color in each spectrum. 
This procedure does not ensure that the calibration is good in the complete wavelength 
range; therefore we multiply by a smooth spline with knots at the effective wavelength of 
the SDSS filters until the synthetic u — g,g — r,r — i and i — z colors match the observed 
colors obtained, again by using polynomial interpolation of the light curves. 



The i^'-corrections for the same filter (IHamuy et al.lll993l ) calculated from the modified 



spectrophotometrically calibrated spectra using SDSS passbands are listed in Table El The 
i^'-corrections are probably accurate to ±0.05 mag for g and r filters, and to ±0.1 mag for 
u, i and z where we had to extrapolate the spectra. However, our estimate is not precise 
and we can not exclude the possibility of even larger errors. After fitting with a low order 
polynomial we use the results to transform the observed SDSS magnitudes in Table [2l to the 
rest-frame. 

We corrected the CSP u'g'r'i' magnitudes to rest-fr ame SDSS uqr i mag nitudes using 



cross-filter i^-corrections according to the prescription of iNugent et al.l (120021 ). These take 
into account the difference between the CSP and SDSS passbands convolved by the SED of 
SN 2005gj and allow us to put all the rest-frame magnitudes in the same system. We find 
that the differences between the same- and cross-filter i^'-corrections are small (< 0.03 mag) 
at all epochs. 

Figure [61 shows the evolution of the rest-frame colors of SN 2005gj as a function of time 



- 13 - 



after explosion. We have corrected the magnitudes by ii'-corrections and Galactic extinction 
in the line of sight. As a comparison we also plot th e color evolution of t he overluminous 
Type la SN 1991T, a typical Type Iln SN 1999el Jpi Carlo et al.l l2002h (both obtained 
from spectral templates of P. Nugent § ), and two previous cases that are thought to be 
Type la explosions in a very dense env ironment: SN 2002i c and SN 1997cy. For SN 2002ic 
we use the published BVI photometry (IHamuy et al.ll2003l ) corrected by Galactic extinction 
{E{B — V) = 0.06; SFD). We calculate cross-filter i^-corrections using th e calibrated spectra 
of SN 2005gj to transform observed magnitudes of the SN at z = 0.0667 (IKotak et al.ll2004l ) 
to rest-frame SDSS magnitudes: B ^ u, V ^ g, I ^ i. The time o : 
is assumed to at ^JD 2,452,581.5 (2002, Nov. 3 UT: lOeng et al 



frame colors in the SDSS system for SN 1997cy (IGermanv et al 



2000 



explosion of SN 2002ic 



20041. We obtain r est- 



Turatto et al.ll2000l) by 



transforming the i^-corrected VRI magnitudes in [Germany et al.l (120001 ) to gri magnitudes 
using S'-corrections calculated directly from spectra of SN 1997cy available online in the 
SUSPECT database We supplement this with synthetic colors from the spectra. We also 
correct by Galactic extinction in the line of sight {E{B — V^) = 0.02; SFD) and assume the 



time of explosion of SN 1997cy to be JD 2,450,582.5 (1997, May. 14 UT; ICermanv et al. 



20001 ). which is very uncertain, since it is taken as the time of detection of the gamma- 
ray burst GRB 970514, which may not have been associated with the SN. Magnitudes for 
SN 2002ic and SN 1997cy are not corrected by extinction in their host galaxies, which is 
unknown. 

Initially the evolution in rest-frame u — g and g — r colors of SN 2005gj, up to ~30- 
40 days after explosion (~10-20 days after maximum for a SN la), is roughly consistent with 
the colors of SN 1991T but ~ 0.2 mag redder in g — r, and evolves to redder colors at later 
times. SN 1991T reaches its maximum colors of {u — g)=1.5 mag and {g — r)=1.0 mag at 
~50 days (30 days after maximum), and after that it enters the nebular phase and becomes 
bluer. At late times the u — g color of SN 2005gj has a slow linear increase and becomes 
systematically redder than SN 1991T at > 70 days, while the g — r color stays approximately 
fiat at [g — r)=0.5 mag (and bluer than SN 1991T) between 60 — 110 days. The evolution 
m u — g color of SN 2002ic is very similar and consistent with SN 2005gj. SN 1997cy has a 
similar color evolutioii but it is ~0.2 mag bluer in g — r color between 60 — 100 days. However, 
Germany et al. J2000h give /^-correction errors of ~0.15 mag and we have further applied 
S'-corrections, therefore this does not imply a significant difference between SN 1997cy and 
the other two Type la/IIn supernovae. 



^http : //supernova. Ibl .gov/~nugent/nugent_templates .html] 
" ^http : //suspect .nhn. ou. edu/'^suspect/ 1 
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The rest-frame r — i and i — z color evolution of SN 2005gj is very different from 
SN 1991T and closely follows the evolution of a SN Iln. The r — i colors of SN 1997cy are 
also consistent with SN 2005gj. From these comparisons it is clear that the colors, a proxy 
for the temperature of the photosphere, of SN 2005gj and two earlier cases of SN la strongly 
interacting with their circumstellar medium are dominated by the radiation coming from the 
ejecta-CSM interaction. 

In Figure [7] we present the ugri light curves of SN 2005gj in absolute magnitudes. We 
also show the light curves of SN 1991T, SN 2002ic and SN 1997cy obtained from the literature 
and corrected to SDSS rest-frame magnitudes as explained above. 

SN 2005gj has peak ugriz absolute magnitudes in the range —20.0 to —20.3 mag, this is 
~ 0.7—1.4 mag brighter than the overluminous Type la SN 1991T. The u light curve is consis- 
tent with a linear decline after peak luminosity at a constant rate of 0.027±0.001 magday"^. 
The g light curve has a ~ 20 day plateau with roughly constant luminosity after maximum 
(20 — 40 days after explosion), then the light curve declines linearly between 40 — 100 days at 
0.018 ± 0.001 magday"^ and continues its linear decay at later times, but with a shallower 
slope of 0.007 ±0.002 magday^^. The r and i light curves have a similar plateau shape be- 
tween 20 — 60 days and a constant linear decay at later times of 0.013±0.001 magday"^. The 
change in slope observed in the (7-band at late times is less clear in rz, but still present. The 
secondary maximum present in ri light curves of SN 1991T and other SN la is completely 
absent in SN 2005gj. 

The light curves of SN 2002ic are fainter than SN 2005gj at all times by 0.3 — 0.6 mag 
depending on the filter (0.3 — 0.8 mag brighter than SN 1991T at peak). The initial decline 
rates from maximum of the u and g light curves of SN 2002ic are intermediate between 
those of SN 1991T and SN 2005gj until around 40 days after ex plosion; after day 50 when 



the ejecta-CSM interaction had become dominant in SN 2002ic (IHamuy et al.ll2003l ). they 
closely resemble the decline rates of SN 2005gj in the same bands. The i band light curve 
of SN 2002ic showed definite evidence for a weak secondary maximum, which is again inter- 
mediate in morphology between the strong secondary maximum observed in SN 1991T, and 
the absence of such a feature in SN 2005gj. The gr light curves of SN 1997cy are consis- 
tent with linear decay of ~ 0.008 magday"^, being ~ 0.4 mag brighter than SN 2005gj at 
150 days. Unfortunately, the light curves of SN 1997cy start at ~ 60 days after explosion 
and there is no information near peak to compare with the luminosities of SN 2005gj and 
SN 2002ic. However, if we extrapolate the gr light curves to the time of peak luminosity we 
find that the absolute magnitudes at maximum of SN 1997cy would be within ~ 0.2 mag of 
SN 2005gj. This has to be taken with caution because of the extrapolation and uncertain 
time of explosion of SN 1997cy. 
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5.2. NIR light curves 

We used the CSP JHKg photometry obtained between 59-166 days after the explosion 
(rest-frame) to construct the absolute magnitude light curves of SN 2005gj. The observed 
magnitudes were corrected by Galactic reddening in the line of sight {Aj = 0.108, = 
0.070, and A^s = 0.044 mag) and i^-corrections. We calculated i^'- corrections for the same 
filter using the spectral templates of P. Nugent for the Type Iln SN 1999el (which are derived 
from black-body curve fits to the photometry), because as we showed in § 15.11 the synthetic 
colors obtained from the spectral templates approximate reasonably well the evolution of 
the redder optical colors of SN 2005gj (i.e., r — i and i — z). The values of the i^-corrections 
are consistent with being constant in this time range: Kj ^ —0.12 mag ,Kh ^ —0.14 mag, 

~ 0.16 mag. 

The JHKs absolute magnitude light curves of SN 2005gj are presented in Figure [81 
For comparison we show the NIR light curves of a normal Type la obtained from syn thetic 



photometry of spectral templates and the Type Iln SN 1999el (jPi Carlo et al.ll2002l ). The 



Type la light curves hav e been shifted in magn itudes to match the mean absolute magnitude 



at maximum of SN la (IKrisciunas et al.ll2004l ). SN 2005gj is 1.7 — 3 mag brighter than a 



normal SN la and SN Iln at 60 days after the explosion and declines in a slower fashion at 
later times. A linear decay is a good fit in all NIR bandpasses at this late times (> 60 days 
after explosion), with decline rates of ~ 0.014 magday"^ in J, ~ 0.013 magday"^ in H, 
and ~ 0.011 magday"^ in Ks- These values are similar to the decline rates in the optical ri 
bands. 

Since there are no template images of the host galaxy obtained in the NIR bands, the 
light curves are preliminary and the analysis has to be taken with caution. 



5.3. Bolometric light curve 

The SDSS and CSP magnitudes were used to produce a quasi-bolometric light curve 
of SN 2005gj covering the optical wavelength range from 3000 — 10,000 K{u — * z). We 
corrected the magnitudes by Galactic extinction in the line of sight and i^'-corrections to 
obtain magnitudes in the rest-frame ugriz filters (see § 15.11 for details). We applied small 
corrections to transform the magnitudes based in the SDSS photometric system into the 
AB system obtained from the SDSS website 0. The AB magnitudes derived in this way 
are transformed directly to bandpass averaged fluxes using the definition of the AB system 



|http : //www. sdss . org/drS/algorlthms/f luxcal .html[ 
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( lOke fc Gunnlll983l ) and they are assigned to the frequencies that correspond to the effective 
wavelengths of the S PSS ugriz filters, calculated from the filters using the definition in 
Fukurita et aP Jl996h : 3567, 4735, 6195, 7510, and 8977 A. We use the trapezoidal rule 
to obtain the integrated flux from u to z, this is between Ai = 3340 A and A2 = 9596 A, 
where the limits of the wavelength coverage are obtained from Ai = Ae//,u — AAu/2 and 
A2 = ^eff,z + AAz/2. We extrapolate linearly the u and z light curves at late times to flll 
in the lack of coverage of the SDSS in z, and CSP in u and z bands, including the MDM 
z-band data at these epochs. The integrated fluxes are converted to luminosity assuming 
a luminosity distance to SN 2005gj of 268.5 Mpc and a spherically symmetric distribution 
of the output energy. We present the integrated, quasi-bolometric luminosities from u to i 
(L(^j^j)) and from u to z {L{u^z)) as a function of time in Table [8l 

In order to estimate bolometric UVOIR luminosities we calculate time dependent holo- 
metric corrections to include the energy output of the SN at wavelengths bluer than w-band 
(A < 3340 A) and redder than 2;-band (A > 9596 A). We flnd that black-body distributions 
with different temperatures are a reasonable approximation of the spectral energy distribu- 
tion of SN 2005gj. We use x^-minimization to flt the optical SED with a two parameter 
black-body function: temperature and a multiplicative scaling factor. The scaling factor is a 
combination of fundamental constants and the square of the angular radius of the spherical 
black-body: (Rbb/d)^, where Rbb is the radius of the black-body and d is the distance to 
the SN. We calculate time-dependent bolometric corrections by integrating the black-body 
distributions in the ultraviolet and NIR/IR regions, and converting the integrated fluxes to 
luminosities as explained above. For the CSP data we also include NIR flux densities derived 
from the reddening and i^-corrected YJHKs magnitudes after fltting the light curves with 
low order polynomials. 

At early phases before peak the bolometric corrections account for ~ 53 — 65% of the 
total integrated luminosity, of which 85 — 93% is from the ultraviolet part of the spectrum 
and only 7 — 15% from the NIR/IR. As the supernova evolves, the ejecta expands and 
shocks the circumstellar gas. The energy emitted in the ultraviolet /blue part of the spec- 
trum declines quickly after maximum light and most of the energy is emitted in the optical 
region, coinciding with the appearance in the spectrum of emission/absorption features of 
the intermediate and iron-peak elements (see § 15.41) . Between 25 — 150 days after explosion 
the bolometric corrections account for ~ 32 — 45% of the total output luminosity, with the 
NIR/IR correction dominating completely over the blue/ultraviolet at > 60 days. 

The bolometric UVOIR luminosities, black-body temperatures and radii derived from 
the fits are presented in Table O The uncertainties in black-body temperatures and radii are 
calculated using the diagonal terms of the covariance matrix obtained from the minimiza- 
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tion. We add a 10% error in the distance to the SN to the error in the black-body radii, which 
comes from the randoni and systematic uncertainties in the value of the Hubble constant 
( iFreedman et al.ll200ll : iRiess et al.ll2005l ). The uncertainties in the bolometric luminosities 
were estimated by propagating errors through the trapezoidal integration of the SED, taking 
into account: uncertainties in the photometry, light curve interpolation and fitting. Galactic 
extinction, ii'-corrections, and distance to the SN. To approximately take into account the 
errors introduced by the bolometric corrections we multiplied these values by \pxi when the 
reduced is greater than 1. 

In Figure [9] we show some examples of black-body fits to the optical SED at different 
epochs. At early times, shortly after explosion, the SED is very well fit by a hot ~13000 K 
black-body. The temperature starts to decrease steadily close to the time of explosion to a 
constant value of ~ 6500 K at 60 days after. A black-body is still a reasonable approximation 
of the SED at later times, but the fits become poorer when emission/absorption features start 
to dominate the spectrum, which is represented in the of the fits (see Table [8]). 

Figure [TU] shows the bolometric light curve of SN 2005gj in the top panel and the evolu- 
tion in temperature and radius from the black-body fits in the lower panels. The early data 
of the bolometric light curve are well fit by an exponential rise in luminosity, Lif) oc e°'^''*. 
The time of maximum bolometric luminosity occurs between 6.6-18.8 days after explosion. 
After maximum, the bolometric light curve is very well approximated by an exponential 
decay in luminosity, linear in the logarithmic scale shown in Figure [TUl L(t) oc 6"°"^^* 
(0.014 magda y"^). This is consistent with th e exponential density distribution of the ejecta 
of Type la SN (IDwarkadas fc Chevalierlll998l ). whereas the distr ibution of ejecta around core - 
collapse supernovae is better approximated by a power- law (IChevalier fc FranssonI |2003| ). 
Extrapolating the pre- and post-maximum fits we get a maximum bolometric luminosity of 
L^"^ = 5.6 X 10^^ ergs~^, which is /-^3 t i mes more luminous th an the Type la SN 1991T 
at maximum light (jContardo et al.l l2000l : IStritzinger et al.l 120061 ). Assuming a bolometric 
correction of 50% at maximum, we find that SN 2005gj was ~1.5 times more luminous than 
SN 2002ic. 



5.4. Optical spectroscopy 

In Figure [TT] we show a comparison of the spectra of SN 2005gj with spectra of SN 2002ic 
and SN 1997cy obtained at similar times after explosion. The spectra of SN 2005gj and 
SN 2002ic are very similar at all times. They are characterized by strong and broad 
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Hydrogen-Balmer lines Ha and H/5 in emissio and a blue continuum at early times that 
becomes redder and increasingly dominated by absorption/emission P-cygni profiles from 
Fe-peak ions (e.g., Fe II, Fe III, Ni III, Si II, S II). 



5.4- 1- Classification 



Benetti et al.l (120061 ) proposed that SN 2002ic-like events are well explained by the core- 
collapse of a massive star in a dense medium, casting doubt in the previous classification 
of SN 2002ic as a Type la supernova. Th e authors found relatively good agreement at all 
times between the spectra of SN 2004aw (ITaubenberger et al.ll2006l ). a Type Ic supernova, 
and SN 2002ic. 



We used the SuperNova I Dentification code, SNID (IMatheson et al. I I2OO5I : iMiknaitis et al 



20071 : iBlondin fc Tonryl 120071 ) . to find the spectra that best match SN 2005gj at different 
epochs. SNID cross-correlates an input spectrum with a library of supernovae spectra. In 
the library we included spectra of 5 normal SN la, two 1991T-like objects, two 1991bg-like 
objects, 4 broad-lined SN Ic (or hypernovae), and 3 normal SN Ic (including SN 2004aw), 
that were chosen to span a wide range of observed properties of SN la and SN Ic. In Table [9] 
we present the supernovae and the epochs of the spectra in the library. We fixed the redshift 
of SN 2005gj at z = 0.0616 and allowed for a range around the mean redshift of Az = 0.02 
to find cross-correlation peaks. The Balmer lines in emission were clipped from the input 
SN 2005gj spectra to avoid spurious cross-correlation signal with library spectra that contain 
emission lines from the host galaxy. 



Figure [T2] shows the spectra of SN 2005gj at four epochs and the library spectra with 
the highest cross-correlation significance from SNID. Type la supernovae spectra are a better 
match to SN 2005gj at most epochs, with 20 of the 26 (77%) epochs having best matched a 
type la spectrum (45%-91T, 45%-normal, 10%-91bg) and a similar fraction for the next-best 
matches. The broad-lined SN Ic 1997ef and 2002ap are the best matching spectra for 6 
epochs, all of them between 26-46 days after explosion. We repeated the sa me procedure 



using five spectra of SN 2002ic obtained between 24-84 days after explosion (IHamuy et al. 



20031 ). All the spectra of SN 2002ic are be tter matched with SN la, in contradiction with 
the results obtained by lBenetti et al.l (120061 ). 



The continuum of SN 2005gj is well approximated at all times by the sum of a scaled 
spectrum of the overluminous Type la SN 1991T, at the same epoch after explosion as 



^As shown in Figure [2l H7 is also visible in the early spectra of SN 2005gj. 
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SN 2005gj, and a fourth order polynomial. A normal SN la does not fit as w ell as SN 1991T. 
This procedure is ver y similar to the fits to SN 2002ic (IHamuy et al.ll2003l ) and SN 2005gj 
( lAldering et al.ll2006l ) presented in previous studies. In Figure [T5] we show examples of the 
spectra decomposition at four epochs. We excluded from the fit a region of ±100 A around 
the Ha and H/3 lines and obtained a good fit for the remainder of the spectrum. 



5.4-2. B aimer lines 

We analyzed the Balmer emission features in the spectra using the sum of two Gaussian 
components to model the line profiles. This decomposition gives r nuch better fits for Ha 
at all epochs than a single Gaussian and it is physically motivated ( Chugai 1997a| jbl). The 



spectra of Type Iln SN show Balmer features with both a narrow and broad component that 
can be explained by radiation coming from different regions of the ejecta/CSM, whether it is 
direct emission from the shock-heated GSM (broad component) or emission from un-shocked 
gas photoionized by the SN radiation (narrow component). The H/3 line is unresolved or 
only marginally resolved for most of the spectra. Therefore a single Gaussian component was 
used to fit the line profile. We used a third order polynomial to model the local continuum 
around each line that was included in the Gaussian fits. It is important to stress that at late 
times there is a broad Fe II feature intrinsic to the supernova spectrum in the region of Ha 
(see spectra in Figure [T5]) that makes the definition of the continuum less reliable and may 
affect the line measurements. 

The results of the Gaussian fits to the Ha and H/3 emission features, integrated fluxes and 
FWHM, are shown in Table [TD] as a function of epoch of the spectra. We have excluded the 
two spectra with better resolution because they show P-Gygni profiles (see below). We used 
the flux calibrated spectra corrected to match the observed g, r magnitudes (as explained 
in § 15.11) and corrected for Galactic reddening in the line of sight. The FWHM of the 
Gaussian profiles listed in Table [10] were corrected by the resolution of the spectrographs 
(from Table E]). We do not present the values when the component is unresolved. The errors 
quoted for the integrated fluxes are obtained by adding in quadrature an estimate of 10% 
error assigned to the absolute flux calibration and the rms deviation of the Gaussian flts. 

Figure [T3] shows the evolution in time of the FWHM (top left panel), Ha and H/3 
luminosities (top right and bottom left panels) and the Balmer decrement (bottom right 
panel). The FWHM of Ha varies between ~ 130 — 500 km s~^ (narrow component) and 
~ 1200—3800 km s~^ (broad component), with the broad component showing a slow increase 
in time. The FWHM of H/3 varies between ~ 470 — 1700 km s~^ and does not show evident 
evolution. The luminosities of Ha-narrow and B.j3 lines evolve in a similar fashion, increasing 
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at early times to peak at ~ 12 days with luminosities 6 — 6.5 x 10^° ergs~^, then they decay 
and stay roughly constant after 50 days. The evolution of Ha-broad is similar during the first 
50 days, peaking at 1.1 x 10^^ ergs~^, but it shows an increase at later times. Compared with 
the Ha luminosities observed in SN 2002ic, both components are ~ 4 times more luminous. 

The Balmer decrement, the ratio between Ha (sum of narrow and broad components) 
and H/3 fluxes, stays approximately constant during the first 30 — 40 days (mean = 2.5 and 
rms = 0.5) and is cons istent with the theoretical value in Case B recombination of Ha/H/3 = 
2.86 JOsterbrocklllQSQi l At later times it shows an steady increase, reaching Ha/H/3 ~ 7 — 13 
at ~ 80 days. In Case B recombination a Balmer decrement Ha/H/? > 2.86 is usually 
interpreted as evidence for the presence of internal extinction in the host; however, the 
large values observed at late times would produce an Na I D in terstellar absorption d oublet 
easily detectable in the spectra, which is not observed (s e e also Aldering et al. (j2006 )). and 
the evolution in time is not expected. lAldering et al.l (120061 ) proposed that the H level 
populations are in Case C recombination, where the optical depth in the Ha line is high 
implying high densities and greater importance of coUisional processes. In this scenario, the 
observed change in the Balmer decrem ent could indicate that coUisional excitat ion becomes 



i ncreasingly impo rtant at later times (IBranch et al.lll981 



( Deng et al.ll2004r) and other SN Iln, like SN 1988Z flAretxaga et al 



Turatto. et al 



19931 ). SN 2002ic 



19991 ) and SN 1995G 



( iPastorello et al.ll2002l ). have also shown large values of the Balmer decrement and therefore 
may have similar physical processes affecting the formation of the Balmer lines. 

In Figure [T5] we show the regions around Ha and H/3 features in the best resolution 
spectra from ISIS and LDSS-3, obtained at 44 and 115 days after explosion, respectively. We 
clearly detect P-Cygni profiles in all these features, which indicates the presence of an outflow 
moving at ~150-200 km s~^; however, these measurements are limited by the resolution of the 
spectra between ~130-180 kms~^ (FWHM). After correcting for the resolution we obtain an 
outflow velocity of 60-70 kms~^. The detection of P-Cygni-like absorption rules out an H II 
region in the line of sight that could be producing the narrow emission/absorption features; 
the line profiles are intrinsic to the SN. lAldering et al.l (120061 ) detected P-Cygni profiles in 
He I, Ha and H/?, in a high resolution spectrum obtained with LRIS+Keck 71 days after 
the explosion. The y derived a wind velocity of ~ 60 kms~^ consistent with our estimate. 
Kotak et al.l (120041 ) also detected a P-Cygni profile in the a spectrum of SN 2002ic obtained 
256 days after explosion. 
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5.4-3. Line identification near maximum 



We used the parameterized resonance-scattering code SYNOW ( Fisher et alll997l : lFisher 



2OOOI ) to identify the hnes in the spectra obtained near maximum hght of SN 2005gj. SYNOW 
is a fast supernova spectrum-synthesis code used for direct (empirical) analysis of supernova 
spectra, mainly to identify the lines, their formation velocities and optical depths. The 
code is based on simple assumptions: spherical symmetry, homologous expansion, a sharp 
photosphere that emits a black-body continuous spectrum, and line formation by resonance- 
scattering, treated in the Sobolev approximation. We have used the latest version of the 
code that includes a Gaussian distribution of optical depths. 



Figure [16] shows the spectrum of SN 2005gj at 17 days after explosion (2 days before 
g maximum) and the best synthetic spectrum obtained with SYNOW. We also show for 
comparison the spectrum of SN 1991T obtained at -3 days w ith respect to the t ime of 
B maximum. The spectra have been locally normalized as in iJeffery et al.l (120061 ). The 
synthetic spectrum has a black-body continuum temperature Tfyj, = 11000 K, photospheric 
velocity Vphot = 10000 km s~^, and excitation temperature Tf>xc = 10000 K. We find a 
reasonably good match with the spectrum of SN 2005gj using the following lines/ mult iplets: 
Fe III A4404 and A5129, Si III A4561, Ni III, S II A5468 and A5633, and Si II A6355. 
These lines are characteristic of the overluminous and spectroscopically peculiar. Type la 
SN 1991T around maximum light with s trong Fe III featur es and weak S II doublet and Si II 
Jjefferv et al.lll992l : iMazzah et aDll995l : iFisher et al.]ll999h . 



The main discrepancy between the SYNOW modeling of SN 2005gj and SN 1991T is 
in the optical depths of the lines. The fit to SN 2005gj needs unphysically small optical 
depths, approximately 1/lOth of the values used for SN 1991T around maximum light. We 
interpret this as an effect of the extra continuu m radiation that is added by the ejecta-CSM 



interaction, which is veiling (IBranch et al.ll2000l ) the supernova lines (e.g.. lHamuy et al.ll2003 



Aldering et al.ll2006[ ). This interpretation is supported by the good agreement obtained from 
fitting the spectra of SN 2005gj using a simple polynomial continuum added to the spectra 
of SN 1991T at the same epochs after explosion (see Figure [13]). 



6. Discussion 



We have presented extensive spectroscopy and optical/NIR photometry of SN 2005gj ob- 
tained by the SDSS-II and CSP supernova groups during the first ~150 days after explosion, 
and also an X-ray observation at 74 days that gives an upper limit on the X-ray luminosity. 
We have shown the remarkable similarity in spectroscopic and photometric properties be- 
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tween SN 2005gj and SN 2002ic, which is thought to be the first clear case of a thermonuclear 
supernova explosion embedded in a dense CSM. The observational properties of SN 2005gj 
support this interpretation, they are summarized as follows: 

• Spectroscopic evidence for a shock propagating into an Hydrogen-rich medium close 
to the site of the explosion inferred from the presence of Balmer lines with narrow 
(FWHM~ 200 - 500 kms"^) and broad (FWHM~ 1500 - 3000 kms"^) components 
at all times. The Balmer lines show P-Cygni profiles in the highest resolution spectra 
obtained at 44 and 115 days after explosion, these detections show the presence of a 
slow (~ 100 kms~^) moving outflow. Both observations support the interpretation of 
the supernova ejecta interacting with a dense circumstellar material. 

• Spectrum evolves from a very blue continuum (13000 K black-body) similar to SN Iln 
at ~7 days after explosion to a redder continuum at later times with P-Cygni absorp- 
tion/emission profiles. The strongest lines present around maximum are identified with 
singly and doubly ionized iron-peak elements (especially strong Fe III, weak S II and 
Si II) and the spectra are well matched by the overluminous Type la SN 1991T diluted 
with a polynomial continuum at similar times after explosion. 

• Very luminous and slowly declining bolo metric light curve. The linear decay in lumi- 
nosity after peak (~ 0.014 magday"^) suggests an exponential density distribution of 
the ejecta, which is consistent with the eject a- density profiles obtained from simulations 
of SN la. 



The data presented here on SN 2005gj makes the interpretation of 2002ic-like event s 
as thermonuclear superno yae in a dens e CSM , initially proposed by iHamuy et al.l (120031 ) . 
stronger. In contrast with lBenetti et al.l (120061 ) . we find that the overall shape of the spectra 
of SN 2005gj are more consistent with spectra of SN la at different epochs. Specifically, 
Type Ic SNe usually do not show the S II doublet at ~ 5500 A around maximum light; in 
fact SN 2004aw sho^ys only two very weak notches at the wavelengths of S II near maximum 
(ITaubenberger et al.ll2006l ). This is on e of the identifying features in SN la spectra, also 
present in the overluminous SN 1991T (jPhillips et al.lll992l ). In the spectrum of SN 2005gj 
obtained at 17 days (see Figure [T6|) we detect a weak double absorption that we identify with 
S II, that is much stronger in the spectrum of SN 2002ic around maximum light. We can 
see on the top of Figure [TT] that the spectrum of SN 2002ic obtained 24 days after explosion 
clearly shows this feature. Other conspicuous features observed in SN 2005gj and SN 2002ic 
around maximum are Fe III and Si II. These features are present in SN 1991T, but Fe III is 
not observed and Si II is generally weaker in SN Ic. 
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SN 2005gj has stronger ejecta-CSM interaction than SN 2002ic. The peak bolometric 
luminosity is ~ 1.5 times brighter and the broad and narrow components of Ha are ~ 4 times 
more luminous in SN 2005gj. The fact that the SN 1991T features are weaker in SN 2005gj 
compared with SN 2002ic at similar epochs is consistent with this interpretation, because the 
supernova features are more diluted by the stronger continuum. The absence of evidence for 
a secondary maximum in SN 2005gj, whereas the i band light curve of SN 2002ic does show 
a hint of such a feature, is likewise consistent with the ejecta-CSM interaction in SN 2005gj 
having been stronger than in SN 2002ic. 



6.1. Structure of the CSM 



The circumstellar interaction of core-collap se supernovae with a c i rcum stellar medium 



has been studied in detail in the literature (see iChevalier fc FranssonI (120031 ) for a review). 
When the fast moving ejecta encounters the approximately stationary CSM, a forward shock 
moving into the CSM (also called circumstellar shock) and a reverse shock develops. The 
fast-moving shockwave implies large post-shock temperatures, therefore radiating energy in 
the X-ray regime. The density distribution of the ejecta and the CSM can be well described 
by power-laws in radius, which leads t o a set of self-s imilar analytical solutions for the 
evolution of the shock radius in time (IChevalierl Il982l ). The physics of the ejecta-CSM 
interaction in the case of thermonuclear supernovae is basically the same, the main difference 
i s in the distribution of the eje cted material which follows an exponential function in velocity 
( iDwarkadas fc Chevalienl 19981 ). In this case the solutions are no longer analytic. The density 
profile of the shocked region is different in the case of exponential ejecta expanding into a 
constant density medium, but the similarity increases for expansion into a wind profile whose 
density decreases as oc r^^. 

A simple self-similar model of a SN shock expanding into a medium with a power- 
law density decline, as suggested for core-collapse SNe by IChevalierl (119821 ). is ruled out 
for this object by several observations: the exponential decrease in luminosity, suggesting 
an exponential ejecta density profile; the strange behavior of the broad and narrow Ha 
components; and the decrease in the blackbody radius at later times. Detailed calculations 
of the SN-CSM interaction would require highly detailed hydrodynamic modeling, which 
are beyond the scope of this paper. Instead herein we focus on trying to explain the basic 
features of SN-CSM interaction as deduced from the observational data. 

The initial velocity of a SN shock wave as it breaks out from the surface is at least of the 
order of 2 x 10^ kms~^. The broad Ha velocities that are seen in the first week or so are of 
the order of 1500 kms~^, and increase to more than twice this value after ~ 50 days. These 
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velocities are almost an order of magnitude smaller than expected SN blast wave velocities 
in the early stages, and a factor of few smaller even after ~ 50 days. Furthermore, the 
SN shock velocity would be expected to gradually decrease as the shock moves outwards, 
whereas the Ha profile actually indicates an increasing velocity after ~ 50 days. 

For these reasons, we suggest that the broad Ha lines do not indicate the SN velocity. 
Instead, we put forward a scenario of a shock expanding into a two-component ambient 
medium: a low density wind in which are embedded high-density clumps. In this picture, 
there should theoretically exist three different velocity components: a broad velocity com- 
ponent, which is not easily seen in this case, and is related to the velocity of the blast wave 
itself; an intermediate velocity component, which is what we have referred to as the broad 
Ha and is related to the velocity of the shock driven into the clumped material; and a 
narrow velocity component, which may be related to the narrow Ha and is representative 
of the velocity of the amb ient medium. This scenario is like the scenario put forward by 



(jChugai fc Danzigerlll994l ) to explain the origin of the broad, intermediate and narrow line 
components in SN 1988Z. The large Ha luminosity of SN 2005gj at late times is very sim- 
ilar to that seen in other Type Iln SNe, and is especially large considering that this was 
a Type la. However, there are significant differences. We do not see a really broad line 
component representative of the SN velocity, although there are some suggestions that this 
may be appearing at late times. In particular, the Ha profile of the spectrum obtained at 
~ 150 days is better fitted by three components, including a very broad component with 
FWHM ^ 7000 kms-\ 

Our scenario envisions the Type la SN shock wave expanding in a clumped medium 
presumably formed by mass-loss from a companion star. The broad component is not easily 
visible in Ha initially because the forward shock is not radiative. The density of the clumps 
is much higher than that of the interclump (ambient) medium. When the SN shock wave 
interacts with a dense cloud or clump, it drives a strong shock in to the clump. A reflected 



shock is driven back into the expanding ejecta (iKlein et al.lll994l ). Assuming pressure equi- 
librium, the ratio of the velocity of the clump shock to that of the blast wave is inversely 
proportional to the square-root of the ratio of the clump density to that of the interclump 
medium. The optical emission arises from behind the clump shock, probably by reprocessing 
of the X-ray emission. 

In this model, the intermediate component represents the velocity of the clump shock, 
which is probably radiative. If we assume that the initial velocity of the SN shock wave 
is ~ 20000 kms~^ and the broad Ha emission velocity is ~ 1500 kms^^, then the ratio of 
velocities is 13 — 14. This indicates that the clump density is about 14^, or ~ 200 times the 
interclump density. Note that the optical emission, which goes as density squared, will then 
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be 200^ times, or about 40,000 times greater compared to that from the interclump medium. 
This is consistent with the fact that no broad hne emission is seen from the interclump 
medium. If the initial velocity is much higher, as is conceivable, the clump density could 
be up to ~ 50% higher, and the ratio between the emission from the dense clumps and 
interclump medium even larger. 

What value of the clump density is suggested? A shock wave traveling at 1500 kms~^ 
would be radiative if it were expanding in a medium whose density is greater than ~ 
10^ cm~^, whereas a 2500 kms~^ shock would require minimum densities of the order of 



10^ cm~^ (jPraine &: McKed Il993l ) in order to be radiative. The CSM density, being two 
orders of magnitude smaller, would then to be > 10^ cm~^. These are just minimum values, 
and it is conceivable t hat the actual c l ump density is much higher. This result is consistent 



with the conclusion of lAldering et al.l (120061 ). 



The observations show that the broad Hct width increases after 50 days, suggesting 
an increase in the clump shock velocity at later times, which could perhaps be due to a 
decrease in the clump density. Conversely, however, the luminosity of the Ha also increases, 
suggesting an increase in the electron density. At the same time, we would expect the SN 
shock to be decreasing in velocity as it continues its outwards expansion. 

We suggest that the way to reconcile these observations is a scenario in which the density 
within the clump medium starts out higher than 10^ cm~^, probably as high as 10^° cm~^ in 
the first few days, and decreases gradually outwards. The almost constant behavior of the 
FWHM of the broad Ha suggests that the density profile of the ambient medium is flatter 
than r~^. Since we want the clump shock to be radiative even when the shock velocity 
is almost 3000 kms~^, this suggest that the density at ~ 150 days is greater than about 
10^ cm~^. And since the density is decreasing outwards, we infer that the density close in is 
even larger. Over the entire period of observations the clump density is large enough that 
the shock driven into these clumps is always radiative. The density of the ambient medium 
is two orders of magnitude smaller, as discussed above. The high bolometric luminosity is 
consistent with these values. 

For the first ~ 50 days the Ha emission arises only from the radiative shock driven into 
the dense clumps. However, by ~ 50 days the SN forward shock, which is decreasing in 
velocity, enters the radiative regime, and the cooling shell of material begins to contribute 
to the Ha luminosity. The velocity of the SN shock is quite large, and its contribution 
initially is not a large fraction of the total Ha luminosity. But as it expands outwards, 
its velocity decreases and the shock becomes more radiative, and the contribution to the 
total Ha luminosity increases, more than compensating for the decreasing density. If this 
conclusion is correct, then we would expect that a broad velocity component would be visible 
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in the Ha spectra, whose intensity would gradually increase with time even as the FWHM 
decreases. Although the underlying supernova contamination makes it hard to isolate a 
broad component, it is suggestive that by day ~ 150 the spectrum is best fit by a third, 
much broader component of the velocity, thus providing support for this line of reasoning. 

Finally, in this model the narrow line emission arises from the unshocked slowly expand- 
ing ambient material, presumably the outflow that we find expanding at ~60 km s~^. We 
note that although the width of the narrow line Ha emission as listed is higher, it is still 
unresolved, and it is possible that within the limits of resolution the narrow line component 
and outflow velocity are indeed the same. 

To summarize, in this model the Type la SN expands in a clumped ambient medium, 
with the clump density about ~ 200 times that of the surrounding medium close in to 
the star, and decreasing as we go outwards. The Ha emission initially arises mainly from 
the shock driven into the dense clumps. The SN shock propagating into the interclump 
medium begins to enter the radiative regime around day 50, and its contribution to the Ha 
emission gradually increases beyond that coming from the clumped medium, leading to the 
gradual rise in the Ha em ission. We note tha t several features of this model are similar 
to the model presented by IChugai et al.l (120041 ) for SN 2002ic, thus further supporting the 
similarity between the two supernovae. 

The upper limit on the X-ray luminosity obtained at 74 days after the explosion can 



put a constrain on the mass loss rate from the precursor or companion (e.g., Ilmmler et al. 
20061 ). Assuming that the X-ray luminosity is dominated by em ission from the r everse shock 
we obtain M < 2 x 10~^ M0yr^^ (2cr) using Equation 3.10 in iFransson et al.l (119961 ). This 
value has to be taken as an approximate estimate because we are making several assumptions 
about the physical properties of the ejecta-CSM interaction that should be calculated using 
detailed hydrodynamical simulations: a constant velocity of the shock, Vgh ~ 8000 kms~^; 
a solar composition of the CSM material; an electron temper ature at the reverse shock of 
Te = 10'' K, which comes from the modeling of SN 2002ic (INomoto et al.ll2005l ): a fiat 
density profile of the CSM, p cx r~^; and a power-law ejecta density profile with index n = 7 
dNomoto et allll984j l 



We can also estimate a mass loss rate from the companion using the density of the 
ambient medium (n ~ 10^ cm~^), the initial optical radius of the CSM {R f» 10^^ cm), and 
the velocity of the wind: M = 4tt R"^ v p, this is assuming a fiat density profile for the CSM. 
We obtain: M ^ 2 x 10^"^ Mq yr^^ which is in agreement with the 2a upper limit calculated 
from the X-ray luminosity. 



The presence of Balmer lines in emission in the first spectrum obtained 6.6 days af- 
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ter explosion shows t hat the ejecta started to interact with the CSM at an earher epoch 
( lAldering et aLll2006l ). Extrapolating linearly to zero flux the early increase of Ha and H/5 
flux;es we find that the ejecta-CSM interaction started 3 ± 1 days after explosion, which 
gives an internal radius of the CSM Ri ^ 1.1 x 10^^ cm. The outer radius of the CSM 
can be estimated assuming a constant velocity of the shock of Vsh ~ 8000 kms~^ over the 
first year. We detect Ha in emission in a spectrum obtained at 368 days after explosion, 
which will be presented elsewhere, putting a lower limit on the outer radius of the CSM, 
i?o > 3 X 10^^ cm. This is also consistent with a Type la SN with an exponential ejecta 
density profile expanding outwards in a medium of average density > 10^ cm~'^. 

In the interpretation above we assume that the broad component of the Hydrogen 
Balmer lines originate in the dense clumps, while the narrow component arises from the 
photoionized un-shocked gas. However, Thompson scattering of the lines has been considered 
as an alterr iative mechanism that can explain relatively well the sy mmetric line profiles of 
SN 2002ic JWang et al.ll2004h and SN 2005gj jAldering et all 120061 ). In this scenario, both 
components would arise from a single high density region. The total mass of hydrogen in 
the emitting region would be f» 2 X lO'^ (lO^V^e) M©, where Ue is the average electron 
density in the emitting zone, as calculated from the luminosity of the Ha line at maximum 
using the Case B recombination coefficient. The electron density must be sufficiently high. 



10^° cm '\ to be consistent with the line ratios of He lines observed in the spectra 



( jAldering et al.ll2006r). and a high e lectron density would explain the non-detection in X-ray 
and radio (ISoderberg &: Fraill l2005l ) . However, it is unlikely that this model would be able 
to explain the initial constancy and then rise of the broad Ha luminosity. 



6.2. Rates, hosts galaxies and possible progenitors of SN 2002ic-like 

supernovae 



The SDSS-II Supernova Survey has a well understood discovery efficiency of SN la at 
low redshift {z < 0.1), which allows us to obtai n an accurate supernova rate measurement 
controlling systematic errors (jPilday et al.l 120071 ). In the fall 2005 season there were a total 
of 16 spectroscopically confirmed SN la (including 1991T-like and 1991bg-like objects) at 
z < 0.12, on e photometric ident ification, and the spectroscopically confirmed peculiar events: 
SN 2005hk Jphilhps et al.ll2007h and SN 2005gj. Since the detection efficiency of 2002ic-like 
objects has not been carefully modeled, we can only put a lower limit on the fraction of these 
events. The spectroscopic confirmation of one object at z < 0.12 puts a lower limit of 5^4% 
(68% confidence) in the fraction of 2002ic-like events among SN la at low redshift. From the 
previously known (2002ic) and probable events (1997cy and 1999E) the estimated fraction 
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IS ~ 



1% of SN la discovered between 1997 and 2002, which is consistent with our hmit. 



In the fall of 2006 we obtained the spectrum of a slowly declining supernova that was 
discovered in 2005, but did not have a spectroscopic classification, SN 7017[3. To our surprise, 
the late spectrum of SN 7017 resembles that of SN 2005gj one year after explosion and the 
early photometry also shows similarities which lead us to classify it as the highest redshift 
SN 2002ic-like object observed to date, ai z = 0.27 (Prieto et al. 2007, in preparation). 
Considering SN 7017 in the spectroscopically confirmed sample of SN la during the 2005 
season, a total of 129 SNe at ^ < 0.42, we have that 2/130 (1.5%) are SN 2002ic like objects, 
which is consistent with the low limit on the fraction at low redshift estimated before. 
However, this fraction has to be taken with extreme care and probably does not reflect 
the true fraction. This is because the discovery efficiency of SN la declines as a function 
of redshift and the total number of spectroscopically confirmed SN la does not include SNe 
with good la-like light curves that were not spectroscopically classified. A more careful study 
of the rates of SN 2005gj-like supernovae in the SDSS-II is planned for a future publication. 

The host galaxies of supernovae can provide important clues about their progenitors. 
The host of SN 2005gj is a very blue, l ow-luminosity dwarf ( Mr ^ —17), and has an irregular 
morphology with no well defined core. lAldering et al.l (120061 ) combined the SDSS photometry 
with UV imaging from GALEX to construct an SED of the galaxy. They constrained the 
metallicity to Z < 0.3 Zq, with a burst of star formation ~ 200 Myr ago. SN 2002ic has a 
late type (Sbc) spiral host with a well defined core. The host of SN 1997cy is also a blue, low- 
lumin osity {My —18.2), and low surface brightness dwarf irregular galaxy (iGermany et al. 



20001 ). GALEX has imaged the positions of SN 1997cy, 2002ic and 2005gj, and their hosts 
galaxies are all detected in the Near-UV (NUV) band. Their absolute magnitudes in the 
NUV (AB magnitudes) are between -16.6 (SN 2005gj) and -17.3 (SN 1997cy and 2002ic). 
They are low-luminosity late type gala xies, ~ 1 — 1.7 mag nitudes fainter than L* galaxies 



observed by GALEX at redshift ^ < 0.1 (IWvder et al.ll2005f ). The host gal axy of SN 1999E i s 
a late spiral with a nuclear starburst first observed by the IRAS satellite (lAUen et al.lll99ll ). 
From 2MASS photometry, its abso lute magnitude is 1 mag brighter than an L* galaxy in 
the K-hemd (IKochanek et al.ll200lf ). SN 7017 at redshift z = 0.27, has a blue, dwarf-like 
host galaxy with absolute magnitude in B of —17.9. 

The host galaxies of the five SN 2002ic-like objects known share some common prop- 
erties: they are late type galaxies, irregulars and late spirals, most likely with recent star 
formation. Four of the host galaxies have low luminosities, similar to the Magellanic clouds. 



®This is the internal name given by the SDSS-II CoUaboration. It was not announced in an lAU circular 
because of the late spectroscopic classification. 
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which indicates they are low metaUicity systems. For example, a dwarf galaxy with intrinsic 
lumin osity Mb = —18 has an Oxygen abundance of 12 + log(0 /H) ^ 8.4 (Ivan Zee et al. 



2006 ). which corresponds to 1/3 the solar Oxygen abundance (jPelahaye &: Pinsonneault 



20061 ). On the other hand, the host galaxy of SN 1999E has a K-hand luminosit y, that when 



conver ted to metaUicity using the luminosity-metallicity relationship derived by lSalzer et al. 



( I2OO5I ). makes it consistent with the solar value. The host luminosities are only an approx- 
imate indicator of their metallicities, therefore spectra of the hosts are needed to infer the 
metallicities and star formation rates (SFRs) of these galaxies. However, it is interesting 
to note that the range of host galaxy properties of SN 2002ic-like e vents seem to be in con- 
sistent with the host galaxi es of GRBs associat ed with supernovae (IStanek et al.ll2006l ) and 
broad-lined type Ib/c SNe JModiaz et al.lboOTh . 



Type la supernovae are observed in all types of galaxies. There is a well established 
correlation between the morphology of their host galaxies and the peak luminosity of the 
SNe: brighter supernovae (1991T-like) tend to explode in late type spirals and irregulars with 
recent star formation, while intrinsically fainter ev ents (1991bg-like) a r e observed mainly in 



early type galaxies with an old stellar populat ion (jHamuy et al.lll995l . Il996at iBranch et al. 



19961 : iHamuy et al.ll2000l : iGallagher et al.ll2005l ). This environmental effect and observations 



of the local supernova e rate as a f unction of host galaxy p r operti es ( ICappellaro et al.lll999 



Mannucci et al.ll2005l ). motivated IScannapieco &: BildstenI (120051 ) to parametrize the delay 
time distribution, time between star formation and the appearance of SNe, and the rates 
with a two-component model having a piece proportional to the SFR of the host galaxy (or 
prompt, they explode ~ few x 10^ yr after an episode of star formation), and a second piece 
proportional to the total stellar mass {delayed component, they explode on scales of a few 
Gyr after the onset of star formation). The difference in age of the stellar populations of these 
subclasses suggests that the progenitors may also be different: prompt SN la would come 
from more massive progenitors. The host galaxies of all five SN 2002ic-like events known 
are broadly consistent with the properties of the hosts of prompt SN la, which suggest a 
real association given that the best studied SN Ila to date, SN 2002ic and SN 2005gj, have 
spectral characteristics similar to 1991T-like events. 

Several p r ogenit ors have been discussed in the literature for SN 2002ic and SN 2005gj. 
Livio &: Riesd (120031 ) proposed that SN 2002ic is a rare case of a double-degenerate binary 
system, a white dwarf (WD) and the core of an AGB star spiraling-in through gravitational 
wave losses, in which the explosion occurs during or immediately after the common-envelope 
phase (a few hundred to a few thousand years of duration). The difference in line strengths 
of the Balmer em ission lines observed for SN 2002ic and SN 2005gj makes this scenario 
unlikely. Also, as lAldering et aD J2006h points out, in both SN 2002ic and SN 2005gj the 
mass loss stopped only a few years before explosion, which is too short compared with the 



- 30 - 



timescale for gravitational wave radiation to produce the merger of the core and the WD. 



Another p ossible progen i tor in itially proposed by iHamuy et al.l (120031 ) and favored by 
the models of IChugai et al.l (120041 ). is the explosion of the Chandrasekhar-ma ss Carbon- 
Oxyg en core of a massive AGB star in a degenerate medium, a supernova Type 1.5 (jiben fc Renzini 
19831 ). where the dense Hydrogen-rich CSM would come from the outer layers of the AGB. 
In order for the core to grow to the Chandrasekhar mass, the radiatively driven winds from 
the AGB have to be weak enough, a condition t hat is only met in a very low-metallicity 
environment like the Galactic halo (jZijlstrall2004l ). At face value, the range of host galaxy 
metallicities for SN 2002ic-like events inferred from the luminosity-metallicity relation does 
not support the SN 1.5 scenario, although admittedly these are average metaUicities and do 
not tell us the actual range of metallicities of the progenitors. 



Han &: Podsiadlowskil (120061 ) proposed that SN 2002ic could be produced through the 
"super-soft channel" , the most common single-degenerate model for the progenitors of SN la. 
In this scenario the white dwarf is accreting material from a main sequence, or slightly 
evolved, relatively massive companion (~ 3 Mq) and experiences a delayed dynamical in- 
stability that leads t o a large amount of m ass-loss from the system in the last few x 10'^ yr 
before the explosion. lAldering et al.l (120061 ) notes that the estimated main-sequence mass of 
the progenitor of SN 2005gj of ~ 2M(7^, calculated using t he age of the starburst of its host 
galaxy, is consistent with the iHan fc Podsiadlowskil (120061 ) model. Also, the predicted frac- 
tion of SN la that would be produced through the "delayed dynamical" channel is 0.1-1%, 
consistent with the limits we have obtain from the detection of SN 2005gj in the SDSS-II 
Survey. 



In general terms, the progenitor model proposed by iHan fc Podsiadlowskil (120061 ) suc- 
cessfully reproduces the observational properties of SN 2002ic and SN 2005gj. However, it 
is still very early in the study of this new sub-class of SN la. It would be interesting to 
see in the near future the results of theoretical modeling exploring other single degenerate 
configurations (e.g., AGB donor) and detailed hydrodynamical modeling of the ejecta-CSM 
interaction of SN 2005gj using the observations of the early photometric and spectroscopic 
evolution presented in this work. 
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Table 1. SDSS ugriz and CSP u'g'r'i' photometry of comparison stars in common in the field of SN 2005gj. 



Star 
ID 


a { 


;j2ooo.o) 


<5 (J2000.0) 




SDSS 


u 


CSP 




SDSS 


9 


CSP 




SDSS 


r 


CSP 




SDSS 


i 


CSP 




z 

SDSS 


1 


03 


01 


09. 


.56 


—00 


33 


52.5 


18 


.639(027) 


18, 


,690(054) 


17, 


,367(018) 


17, 


,362(009) 


16, 


,837(022) 


16, 


,823(013) 


16 


,643(018) 


16, 


,586(010) 


16, 


,549(020) 


2 


03 


01 


06, 


,29 


-00 


32 


59.0 










19, 


,122(027) 


19, 


,043(018) 


17, 


,609(014) 


17, 


,575(015) 


16 


,191(017) 


16, 


.051(017) 


15, 


,403(018) 


3 


03 


01 


14, 


,03 


-00 


31 


48.7 


19 


.843(039) 


19, 


,899(097) 


17, 


.676(016) 


17, 


.653(023) 


16, 


,745(012) 


16, 


,726(015) 


16, 


.404(015) 


16, 


,362(010) 


16, 


.215(017) 


4 


03 


01 


14, 


,44 


-00 


34 


44.3 










20, 


.116(026) 


20, 


.038(045) 


19, 


,579(018) 


19, 


,556(062) 


19 


.356(025) 


19 


.113(067) 


19 


.207(040) 


5 


03 


01 


12, 


,99 


-00 


31 


22.5 






20, 


,262(436) 


18 


.856(021) 


18 


.770(015) 


17, 


,345(009) 


17, 


,323(009) 


16 


.595(014) 


16 


.505(010) 


16 


.145(064) 


6 


03 


01 


21, 


,27 


-00 


32 


44.3 


16 


.074(030) 


16, 


,102(011) 


14, 


.556(010) 


14, 


.544(011) 


13, 


,965(016) 


13, 


,961(015) 


13 


.774(007) 


13, 


.728(015) 


13, 


.711(023) 


7 


03 


01 


07, 


,49 


-00 


30 


51.5 


19 


.915(041) 


20, 


,347(285) 


18, 


.385(019) 


18, 


.384(017) 


17, 


,741(020) 


17, 


,735(010) 


17 


.502(016) 


17, 


.451(013) 


17 


.372(021) 


8 


03 


01 


21, 


,46 


-00 


34 


25.2 


18 


.486(027) 


18, 


,533(030) 


17, 


.214(011) 


17, 


.204(009) 


16, 


,694(017) 


16, 


,685(009) 


16 


.508(018) 


16, 


.473(009) 


16, 


.418(016) 


9 


03 


01 


07, 


,16 


-00 


35 


52.2 


18 


.847(021) 


18, 


,923(041) 


17, 


.458(013) 


17, 


.453(009) 


16, 


,940(011) 


16, 


,936(011) 


16 


.765(018) 


16, 


.732(011) 


16, 


.679(024) 


10 


03 


01 


12, 


,26 


-00 


36 


14.5 


18, 


.285(025) 


18, 


,306(026) 


16, 


.956(018) 


16, 


.935(011) 


16, 


,412(011) 


16, 


,404(009) 


16 


.208(017) 


16, 


.181(009) 


16, 


.125(021) 


11 


03 


01 


02, 


,16 


-00 


35 


21.1 










20 


.227(033) 


20 


.118(058) 


18, 


,717(014) 


18, 


,715(017) 


17, 


.861(021) 


17 


.787(014) 


17 


.347(020) 


12 


03 


01 


24, 


,01 


-00 


34 


27.7 










18, 


.927(020) 


18, 


.860(015) 


17, 


,433(014) 


17, 


,404(020) 


16 


.669(018) 


16, 


.611(010) 


16, 


.243(024) 


13 


03 


01 


21, 


,90 


-00 


35 


57.2 


18 


.510(035) 


18, 


,529(031) 


16, 


.922(017) 


16, 


.897(009) 


16, 


,249(016) 


16, 


,226(011) 


15 


.992(017) 


15, 


.958(009) 


15, 


.868(016) 


14 


03 


01 


01, 


,95 


-00 


36 


08.2 


20 


.030(065) 


20, 


,081(196) 


17, 


.833(019) 


17, 


.819(010) 


16, 


,882(007) 


16, 


,864(013) 


16 


.515(017) 


16, 


.470(009) 


16, 


.306(015) 


15 


03 


01 


26, 


,63 


-00 


31 


51.6 


18 


.748(029) 


18, 


,314(545) 


16, 


,915(010) 


16, 


.888(009) 


16, 


,142(013) 


16, 


,117(009) 


15, 


.845(022) 


15, 


.798(009) 


15, 


.708(015) j 


16 


03 


01 


22, 


,53 


-00 


30 


11.5 


18, 


.767(032) 


18, 


,834(109) 


17, 


.458(013) 


17, 


.454(011) 


16, 


,896(015) 


16, 


,882(009) 


16 


.692(013) 


16, 


.652(009) 


16, 


.599(018) 



Note. — Uncertainties given in parentheses are in thousandths of a magnitude. For the CSP photometry with the Swope they correspond to the rms of the magnitudes 
obtained on four photometric nights, with a minimum uncertainty of 0.015 mag for an individual measurement. 
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Table 2. SDSS and MDM ugriz photometry of SN 2005gj 



JD Epoch* 



-2,453,000 


(days) 


u 


9 


r 


i 


z 


Source 


616, 


,94 






9c: 1 Q_i_n Qfi 


9Q Qi -i_n fii 

ZO. o±IlIU.U± 


Z0.0U^U.04 


99 1 K.'x 

ZZ, 103IU. 00 


01900 


626, 


,91 




91 1 9-i-n 7n 


97 cQ-i-n fifi 
z 1 .oomu.ou 


ZZ.00I1IU.40 


9A Q9^n Q*^ 


99 9Q^1 77 
ZZ. zo^ 1 . i ( 


OJ900 


628, 


,90 




9Q Q0-|_1 AK 


OK 794-1 9Q 

ZO. ( Z3I-L.ZO 


99 Q9-i-n A^ 


zo.yo^u.4o 


91 79^0 K.A 


OJ900 


635, 


,91 




99 nO-t-fl AQ 


9fi "^R-t-l 1 1 


9Q QO-i-n ao 
zo.oz^u.oz 


zo.z / mu.oy 


99 1 ^4-0 K.'i 
ZZ. lomu. 00 


01900 


639, 


,95 


1 Q 


lo.oo 1 ^u^o J 


1 S fiQ9f'n9'^'^ 
J.o.uyz^^uzu ) 


1 Q «01 /"ni Vi\ 
lo.OZl^UlU ^ 


1 S 71 S^ni fi^ 




Oi-'OO 


641, 


,95 


o.o 




1 Q 1 /'ni 1 ^ 


10,Z()01UZU J 


1 9, '\AV^(Vi^ (\\ 


1 Q c;i 1 (r^A^ \ 


01900 


644, 


,89 


u.o 


10191 f'n9'^^ 

iO.J-Z±lUZO J 


17 7Q^^ni I"! 


1 7 Q7Q('n99'l 


1 7 Q'^^/'ni 
1 / .yooi^uj-U j 


1 S 1 A1 ^099"! 


'JDS'? 

01900 


656, 


,94 


17.9 




17.355(026) 


17.255(029) 


17.360(064) 




MDM 


657, 


,90 


18.8 


17.949(039) 


17.343(012) 


17.214(012) 


17.275(016) 


17.438(024) 


SDSS 


663, 


,89 


24.5 


18.149(039) 


17.410(050) 


17.158(012) 


17.163(012) 


17.280(027) 


SDSS 


663, 


,92 


24.5 




17.368(025) 


17.160(020) 


17.248(055) 




MDM 


665, 


,91 


26.4 






17.146(028) 


17.215(054) 




MDM 


666, 


,96 


27.3 


18.282(043) 


17.408(012) 


17.121(012) 


17.133(013) 


17.268(015) 


SDSS 


668, 


,87 


29.2 




17.393(014) 


17.125(026) 


17.186(045) 




MDM 


668, 


,88 


29.2 


18.355(028) 


17.401(009) 


17.127(012) 


17.096(016) 


17.240(016) 


SDSS 


669, 


,97 


30.2 


18.347(045) 


17.430(011) 


17.101(014) 


17.105(031) 


17.231(017) 


SDSS 


670, 


,88 


31.0 


18.386(024) 


17.419(007) 


17.125(013) 


17.092(014) 


17.249(021) 


SDSS 


673, 


,85 


33.8 


18.479(028) 


17.450(009) 


17.123(012) 


17.081(014) 


17.229(015) 


SDSS 


675, 


,85 


35.7 


18.551(025) 


17.477(010) 


17.134(029) 


17.088(019) 


17.182(019) 


SDSS 


676, 


,90 


36.7 




17.464(043) 


17.116(026) 


17.138(055) 




MDM 


680, 


,86 


40.4 


18.663(027) 


17.578(009) 


17.133(011) 


17.054(009) 


17.156(016) 


SDSS 


683, 


,92 


43.3 


18.759(062) 


17.656(034) 


17.220(055) 


17.080(015) 


17.061(066) 


SDSS 


686, 


,86 


46.1 


18.887(047) 


17.695(014) 


17.160(011) 


17.060(010) 


17.138(021) 


SDSS 


687, 


.91 


47.1 


18.937(105) 


17.687(013) 


17.206(024) 


17.042(018) 


17.119(026) 


SDSS 


693, 


,86 


52.7 


19.018(081) 


17.780(009) 


17.207(015) 


17.084(014) 


17.130(018) 


SDSS 


699, 


.85 


58.3 




17.864(060) 


17.271(034) 


17.202(057) 


17.166(024) 


MDM 


699, 


,88 


58.4 


19.083(042) 


17.909(013) 


17.273(009) 


17.124(010) 


17.135(020) 


SDSS 


727, 


,70 


84.6 




18.368(025) 


17.639(043) 


17.551(053) 


17.407(015) 


MDM 


737, 


,75 


94.0 




18.525(038) 


17.791(033) 


17.709(063) 


17.470(057) 


MDM 


739, 


.64 


95.8 




18.540(037) 


17.822(035) 


17.744(065) 


17.548(027) 


MDM 



Note. — Uncertainties given in parentheses in thousandths of a magnitude. 
'^Rest-frame days since the time of explosion (JD 2,453,637.93). 
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Table 3. CSP u'g'r'i' photometry of SN 2005gj 



JD 

-2,453,000 


Epoch 
(days) 




u' 




9' 




r' 




i' 


698.69 


57.2 


19 


159(042) 


17 


884(017) 


17 


243(017) 


17 


106(017) 


699.67 


58.2 


19 


196(041) 


17 


910(017) 


17 


261(017) 


17 


100(017) 


702.68 


61.0 


19 


286(042) 


17 


959(017) 


17 


281(017) 


17 


130(017) 


706.69 


64.8 


19 


312(040) 


18 


031(017) 


17 


345(017) 


17 


173(017) 


712.64 


70.4 


19 


619(060) 


18 


157(017) 


17 


418(017) 


17 


249(017) 


720.66 


77.9 


19 


721(131) 


18 


283(023) 


17 


542(017) 


17 


363(017) 


725.65 


82.6 


19 


964(079) 


18 


365(017) 


17 


602(017) 


17 


412(017) 


728.71 


85.5 


19 


938(104) 


18 


430(017) 


17 


658(017) 


17 


453(017) 


736.62 


93.0 


20 


104(084) 


18 


554(017) 


17 


794(017) 


17 


586(017) 


740.64 


96.8 


20 


229(090) 


18 


606(017) 


17 


856(017) 


17 


641(017) 


741.59 


97.6 


20 


238(103) 


18 


635(017) 


17 


872(017) 


17 


651(017) 


746.61 


102.4 


20 


532(220) 


18 


641(025) 


17 


926(017) 


17 


725(017) 


754.58 


109.9 


20 


509(143) 


18 


770(017) 


18 


021(017) 


17 


803(017) 


761.64 


116.5 






18 


847(017) 


18 


100(017) 


17 


910(018) 


763.56 


118.3 






18 


841(017) 


18 


107(017) 


17 


960(017) 


764.58 


119.3 






18 


859(017) 


18 


143(017) 


17 


947(017) 


768.60 


123.1 






18 


882(019) 


18 


200(017) 


17 


978(021) 


773.55 


127.7 






18 


884(026) 


18 


239(019) 


18 


050(021) 


774.56 


128.7 






18 


946(031) 


18 


221(017) 


18 


079(022) 


783.55 


137.2 






19 


017(018) 


18 


322(017) 


18 


117(021) 


786.53 


140.0 






19 


015(017) 


18 


363(017) 


18 


178(020) 


795.55 


148.5 






19 


099(019) 


18 


416(021) 


18 


227(031) 



Note. — 



Uncertainties given in parentheses in thousandths of a magnitude. 
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Table 4. CSP YJHKs photometry of SN 2005gj 



JD 

-2, 453, 000 


Epoch 
(days) 


Y 


J 




H 




Instrument 


700.71 


59.1 


16.565(015) 


16.484(034) 


16 


253(030) 




Retrocam 


704.68 


62.9 


16.591(015) 


16.537(020) 


16 


315(033) 




Retrocam 


709.66 


67.6 


16.628(015) 


16.550(023) 


16 


271(029) 




Retrocam 


714.58 


72.2 


16.673(015) 


16.594(020) 


16 


389(028) 




Retrocam 


718.65 


76.0 


16.725(016) 


16.658(022) 


16 


364(037) 




Retrocam 


722.62 


79.8 


16.832(016) 


16.725(032) 


16 


490(017) 


16.384(096) 


WIRC 


724.68 


81.7 


16.781(016) 


16.716(028) 








Retrocam 


727.69 


84.6 


16.872(024) 


16.757(036) 








Retrocam 


732.71 


89.3 


16.920(019) 


16.839(025) 








Retrocam 


750.61 


106.1 


17.307(016) 


17.152(016) 


16 


812(016) 




WIRC 


7.5.5. .57 


110.8 


17.337(024) 


17.263(024) 


16 


891(025) 


16.745(039) 


PANIC 


756.59 


111.8 


17.428(016) 


17.269(016) 


16 


928(017) 




WIRC 


757.60 


112.7 


17.423(024) 


17.270(024) 


16 


938(025) 


16.807(042) 


PANIC 


773.57 


127.8 


17.677(016) 


17.476(016) 








WIRC 


776.54 


130.6 


17.538(026) 


17.399(070) 


17 


081(068) 




Retrocam 


777.55 


131.5 


17.648(038) 


17.417(056) 


17 


129(119) 




Retrocam 


782.55 


136.2 


17.513(039) 


17.436(073) 


17 


072(114) 




Retrocam 


783.55 


137.2 


17.766(024) 


17.586(024) 


17 


271(025) 


17.045(042) 


PANIC 


785.53 


139.0 


17.765(025) 


17.570(022) 


17 


192(049) 


17.152(177) 


WIRC 


788.54 


141.9 


17.869(016) 


17.594(016) 


17 


335(019) 




WIRC 


797.52 


150.3 


17.855(061) 


17.563(118) 


17 


313(158) 




Retrocam 


800.50 


153.1 


17.840(041) 










Retrocam 


808.52 


160.7 


17.871(058) 


17.551(079) 








Retrocam 


814.49 


166.3 


17.864(084) 










Retrocam 



Note. — 



Uncertainties given in parentheses in thousandths of a magnitude. 



-42 - 



Table 5. Light-curve parameters for SN 2005gj 



Parameter 


Value 


Time of explosion'' 


637.93 ± 2.02 


Time of Umax 


651.77 ± 0.48 


Time of Qmax 


657.80 ± 1.28 


Time of r-max 


672.49 ± 1.39 


Time of imax 


684.51 ± 1.00 


^max 


17.85 ± 0.05 


Qmax 


17.35 ± 0.01 


^max 


17.12 ± 0.01 


^max 


17.05 ± 0.01 


U° ^ 


17.11 ± 0.12 


if max 


16.94 ± 0.09 


r° 

' max 


16.83 ± 0.07 


i° 

"max 


16.89 ± 0.11 


g ,max 


-20.21 


E{B-V)Gal'' 


0.121 ± 0.019 


A, (Gal) 


0.45 ± 0.07 



»JD-2,453,000 

''Magnitudes at maximum in the 
rest-frame, they have been cor- 
rected by Galactic extinction and 

X-corrcctions. Wc assume a negligi- 
ble extinction in the host galaxy. 



'^Prom Schlegel et al. (1998) 
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Table 6. Spectroscopic observations of SN 2005gj 



JD 


Epoch 




Wavelength 


Resolution'' 


Exposure 


—2, 453, 000 


(days) 


Instrument 


Range (A) 


(A) 


(s) 


644.92 


6.6 




3850 - 


7308 


15 


1200 


646.95 


8.5 


ARC-DTS 


3824 - 


10192 


7 


1800 


650.84 


12.2 


ARC-DIS 


3600 - 


9597 


7 


1000 


655.87 


16.9 


MDM-CCDS 


3823 - 


7283 


15 


1800 


665.92 


26.4 


MDM-CCDS 


3883 - 


7341 


15 


2700 


668.83 


29.1 


MDM-CCDS 


3886 - 


7346 


15 


2700 


676.79 


36.6 


MDM-CCDS 


3882 - 


7338 


15 


3600 


684.73 


44.1 


WHT-ISIS 


3924 - 


8901 


4 


1800 


686.79 


46.0 


MDM-CCDS 


3858 - 


7315 


15 


2700 


698.67 


57.2 


duPont-ModSpcc 


3780 - 


7290 


7 


2700 


699.67 


58.2 


duPont-ModSpcc 


3780 - 


7290 


7 


2700 


700.76 


59.2 


MDM-CCDS 


3933 - 


7391 


15 


2700 


702.73 


61.1 


MDM-CCDS 


3856 - 


7310 


15 


2700 


712.73 


70.5 


MDM-CCDS 


3831 - 


7286 


15 


2700 


722.71 


79.9 


NTT-EMMI 


4000 - 


10200 


9 


2700 


724.66 


81.7 


duPont-WFCCD 


3800 - 


9235 


6 


2700 


725.65 


82.6 


duPont-WFCCD 


3800 - 


9235 


6 


2700 


726.66 


83.6 


duPont-WFCCD 


3800 - 


9235 


6 


3600 


727.67 


84.5 


duPont-WFCCD 


3800 - 


9235 


6 


3600 


728.67 


85.5 


duPont-WFCCD 


3800 - 


9235 


6 


3600 


729.67 


86.4 


MDM-CCDS 


3915 - 


7373 


15 


2700 


737.70 


94.0 


MDM-CCDS 


3909 - 


7368 


15 


2700 


751.60 


107.1 


NTT-EMMI 


3200 - 


10200 


9 


2700 


755.62 


110.9 


MDM-CCDS 


3844 - 


7299 


15 


3600 


759.61 


114.6 


Magellan-LDSS-3 


3788 - 


9980 


3 


3600 


799.52 


152.2 


duPont-WFCCD 


3800 - 


9235 


6 


1200 



Note. — Most of the spectra are the combination of multiple observation, the total 

exposure is given. 

^Average resolution obtained from the FWHM of arc-lamp lines. 
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Table 7. ^-corrections of SN 2005gj 



Epoch 
(days) 


Ku 


Kg 


Kr 






6.6 


0.039 


-0.077 


-0.053 


-0.141 


-0.133 


12.2 


0.162 


-0.097 


-0.005 


-0.196 


-0.047 


16.9 


0.226 


-0.052 


-0.055 


-0.135 


-0.223 


26.4 


0.276 


-0.015 


-0.058 


-0.116 


-0.228 


29.1 


0.255 


0.008 


-0.056 


-0.104 


-0.239 


36.6 


0.296 


0.034 


-0.030 


-0.118 


-0.096 


46.0 


0.251 


0.063 


-0.025 


-0.079 


-0.201 


57.2 


0.195 


0.106 


0.049 


-0.123 


-0.152 


58.2 


0.186 


0.101 


0.013 


-0.068 


-0.216 


59.2 


0.313 


0.090 


0.006 


-0.072 


-0.215 


61.1 


0.241 


0.091 


0.011 


-0.096 


-0.136 


70.5 


0.505 


0.086 


0.055 


-0.113 


-0.054 


79.9 


0.439 


0.112 


0.046 


-0.090 


-0.115 


81.7 


0.328 


0.118 


0.063 


-0.096 


-0.057 


82.6 


0.340 


0.115 


0.063 


-0.104 


-0.102 


83.6 


0.378 


0.105 


0.079 


-0.099 


-0.101 


84.5 


0.326 


0.104 


0.071 


-0.104 


-0.026 


85.5 


0.391 


0.112 


0.078 


-0.110 


-0.036 


86.4 


0.330 


0.124 


0.066 


-0.094 


-0.071 


94.0 


0.249 


0.126 


0.094 


-0.081 


-0.092 


107.1 


0.266 


0.141 


0.082 


-0.084 


-0.089 


110.9 


0.337 


0.138 


0.097 


-0.091 


-0.007 
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Table 8. Derived integrated luminosity and black-body fits. 



Epoch 


logL(u^i) 


logL(„^j) 


log Lboi 




Rbb° 




(days) 


(erg s~i) 


(erg s~^) 


(erg s~i) 


(K) 


(IQlS cm) 




1.9 


42.822 


42.862 


43.285(0.109) 


1 0000/1 1 o*i\ 

12898(1193) 


1 00 /O 1 

1.02(0.16) 


0.7 


3.8 


42.990 


43.029 


43.485(0.105] 


13552(1316] 


1.16(0.18] 


0.3 


6.6 


43.141 


43.180 


43.637(0.098) 


13569(1288] 


"1 00/0 0"l\ 

1.38(0.21] 


0.2 


18.8 


43.344 


43.389 


AO 'i-x o 1 r\ r\r\r'\ 

43.713(^0.096] 


10955(831] 


2.32(0.34] 


1.2 


24.5 


43.335 


43.387 


AO a A A 1 r\ iocr\ 

4o.D44(u.iz5; 


9468(612] 


oT/n A o\ 

2.87(0.42] 


1.2 


27.3 


43.335 


43.388 


43.628(0.135] 


nr\£}(\/ Tort\ 

9069(720] 


3.06(0.50] 


2.2 


29.2 


43.335 


43.389 


43.617(0.140] 


OTTO / COO\ 

8772(683] 


oo/n co\ 

3.22(0.53] 


2.3 


30.2 


43.334 


43.388 


43.617(0.137] 


OTi^"l //^00\ 

8761(692) 


3.23(0.54] 


2.3 


31.0 


43.332 


43.386 


AO c-\ '\ 1 r\ ioo\ 

4o.oii(U.izij 


8697(705] 


3.25(0.55] 


2.6 


33.8 


43.325 


43.381 


43.595(0.140] 


8392(685] 


AO 1 r\ crriN 

3.43(0.59] 


2.9 


35.7 


43.316 


43.374 


43.581(0.160] 


8100(582] 


3.61(0.59] 


2.4 


40.4 


43.299 


43.361 


43.559(0.152] 


T/^o A / crot:;\ 

7684(526) 


no /n f A\ 

3.92(0.64] 


2.5 


43.3 


43.270 


43.338 


AO ^ OC\ { C\ 1 >(0\ 

43.532(0.143) 


TO A /oor'\ 

7284(335) 


A J* / Pro\ 

4.24(0.59) 


0.6 


46.1 


43.269 


43.336 


43.526(0.101) 


7142(477) 


4.36(0.73) 


2.8 


47.1 


43.263 


43.332 


43.523(0.139] 


7115(450] 


4.38(0.70] 


2.1 


52.7 


43.247 


43.317 


43.506(0.123] 


6889(425] 


4.58(0.74] 


2.4 


57.2 


43.219 




AO Aor\lr\ r\c\T:\ 

43.480(0.095] 


6958(247] 


4.32(0.54] 


2.0 


58.2 


43.212 


43.288 


AO A '~70 r\ r\r\r'\ 

43.473(0.096] 


6913(248] 


A A i r\ fcrN 

4.34(0.55] 


2.0 


58.4 


43.215 


43.290 


AO Aor\^r\ ioo\ 

43.480(0.133) 


6656(322] 


yi TO /n Tn\ 

4.78(0.70) 


1.6 


61.0 


43.198 


43.274 


43.458(0.107] 


6840(261] 


4.34(0.56] 


2.4 


64.8 


43.176 


43.254 


43.438(0.087] 


6830(249] 


4.27(0.54) 


2.2 


70.4 


43.133 


43.214 


AO 0(~\ A / r\ I nr»\ 


6570(278] 


A 00 / r\ co\ 

4.33(0.58] 


3.4 


77.9 


43.088 


43.171 


43.355(0.098] 


6633(280] 


4.10(0.54] 


2.8 


82.6 


43.057 


43.143 


AO ooi/r» iTvX 

4o.ozl(u.iz7; 


6380(293] 


A 10/n cro\ 

4.18(0.58] 


4.2 


85.5 


43.038 


43.125 


43.305(0.096] 


6459(282] 


A f\ A / r\ n A\ 

4.04(0.54] 


3.5 


93.0 


42.987 


43.076 


43.256(0.103) 


6360(272) 


3.92(0.53) 


3.6 


96.8 


42.963 


43.054 


43.232(0.101) 


6305(279) 


3.87(0.53) 


3.9 


97.6 


42.955 


43.047 


43.226(0.105) 


6318(279) 


3.83(0.52) 


3.8 


102.4 


42.935 


43.026 


43.205(0.133) 


6372(332) 


3.68(0.52) 


4.4 


109.9 


42.896 


42.989 


43.165(0.117) 


6339(313) 


3.53(0.50) 


4.4 


116.5 


42.859 


42.953 


43.135(0.123) 


6549(350) 


3.25(0.46) 


3.8 


118.3 


42.855 


42.947 


43.128(0.138) 


6566(362) 


3.21(0.46) 


4.0 


119.3 


42.847 


42.940 


43.122(0.117) 


6548(357) 


3.20(0.46) 


4.0 


120.2 


42.839 


42.932 


43.114(0.120) 


6536(352) 


3.18(0.45) 


3.9 


123.1 


42.831 


42.924 


43.104(0.127) 


6533(370) 


3.14(0.45) 


4.3 


127.7 


42.818 


42.909 


43.087(0.155) 


6523(399) 


3.06(0.46) 


5.0 


128.7 


42.809 


42.901 


43.078(0.153) 


6486(399) 


3.07(0.46) 


5.2 


137.2 


42.776 


42.869 


43.043(0.169) 


6404(418) 


2.99(0.46) 


6.3 


140.0 


42.766 


42.857 


43.030(0.197) 


6377(421) 


2.96(0.46) 


6.5 


148.5 


42.738 


42.830 


42.999(0.205) 


6246(438) 


2.94(0.48) 


8.0 



''Black-body temperature from the fits to the broadband photometry; la uncertainty are 
given in parenthesis. 

''Black-body radius from the fits to the broadband photometry; la uncertainty are given 
in parenthesis and include a 10% uncertainty in the distance to SN 2005gj. 

'^X^ per degree of freedom of the black-body fits. 
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Table 9. Library of spectra used in SNID 



SN Name 


Class 


Epochs 


Reference 


1990N 


la normal 


-14, -13, -8, -7, -6, 0, 4, 8, 15, 18, 39 


1 


1991T 


la 91T 


-9, -8, -7, -6, -5, -4, -2, -1, 9, 10, 11, 12, 15, 16, 17, 18, 19, 20, 21, 22, 23, 27, 43, 44, 47, 48, 51, 69, 77 


2, 3 


1991bg 


la 91bg 


1, 3, 16, 18, 25, 32, 33, 46, 54, 85 


4 


1992A 


la normal 


-5, -1, 3, 5, 6, 7, 9, 11, 16, 17, 24, 28 


5 


19941 


Ic normal 


-6, -4, -3, 0, 1, 2, 21, 22, 23, 24, 26, 30, 36, 38 


6 


1997ef 


Ic broad 


-14, -12, -11, -10, -9, -6, -5, -4, 7, 13, 14, 16, 17, 19, 22, 24, 27, 41, 45, 47, 49, 75, 80, 81 


7 


1998aq 


la normal 


-9, -8, -3, 0, 1, 2, 3, 4, 5, 6, 7, 19, 21, 24, 31, 32, 36, 51, 55, 58, 60, 63, 66, 79, 82, 91, 211, 231, 241 


8 


1998bu 


la normal 


-3, -2, -1, 9, 10, 11, 12, 13, 14, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 57 


9 


1998bw 


Ic broad 


-9, -7, -6, -3, -2, -1, 1, 3, 4, 6, 9, 11, 12, 13, 19, 22, 29, 45, 64, 125, 200, 337, 376 


10 


1999aa 


la 91T 


-11, -7, -3, -1, 5, 6, 14, 19, 25, 28, 33, 49, 47, 51 


11 


1999ee 


la normal 


-9, -7, -2, 0, 3, 8, 10, 12, 17, 20, 23, 28, 33, 42 


12 


1999ex 


Ic normal 


-1, 4, 13 


12 


1999by 


la 91bg 


-4, -3, -2, -1, 2, 3, 5, 6, 7, 8, 10, 11, 25, 29, 31, 33, 42 


13 


2002ap 


Ic broad 


-5, -4, 3, 8, 10, 17, 19 


14 


2004aw 


Ic normal 


1, 5, 6, 8, 15, 21, 22, 26, 28, 39, 35, 44, 49, 63, 64, 236, 260, 413 


15 


2006aj 


Ic broad 


-6, -5, -4, -3, -2, -1, 0, 2, 3 


16 



DlLeibun dgut ct al.l ill 9911): (2')|jeffery et al.llll992h: (3')ISchmidt et al.l l ll994): ('4) iLeibundgut et al.l(ll993h: f5')lKirshner et alj 
(6) iMillard et alj l ll999.'): (7) llwarnoto et al.l ipOOOl): (8') iBranch et al.l ||2003|') : f 9) Ijha et all lll999h : (1 0') IPatat et al.l boOlh: fil l 



Ref erences 

lll993^, _ - : ; . -- --- ^ . - . . - 

iGaravini et al.1 (120041') :~( 12) iHamuv et al.1 1I2OO2I' ): fl3'l TGarnavich et al.l ll2004h ~ (li) Ical-Yam et al.l bOO'j) : (15) iTaubenbereer et al.l hood) 
(le't lModiaz et al.1 bom . 



Table 10. Results of the Gaussian fits to Ha and H/3 features 









Ha (narrow) 




Ha (broad) 




H/3 


JD 


Epoch 


FWHM^' 


flux'' 


FWHM^ 


flux'' 


FWHM^ 


flux'' 


-2,453,000 


(days) 


(km s~i) 


(10 erg s cm ) 


(km s"-*-) 


(10~^^ erg s-^ cm-^) 


(km s~^) 


(10~^^ erg s-^ cm-2) 


fiAA Q9 


O.D 






10 ( 


U.Ool^U.UO J 


77R 


U.OU^^U.UO J 




o.O 






14C51 


U.OOl^U.Uo J 


-LoU 1 


A CO ('A Afi^ 


DOU.C54 






u.oyi^u.u / y 




1 oc^/fi 1 '^^ 

L.ZDyyj. io ) 


1 AdO 
14DZ 


n 7C;/'A 1 A\ 

U. / Ol^U. iU ) 


DOO.O ( 


1 Q 
io.y 




U.O / \^U.UD ) 


iOOO 


1 ii/'nii\ 

1 . i i I^U. il ) 


iooy 


u. ( oi^u.uy j 




ZD. 4 




U.OOl^U.UO ) 


iooy 


i.Uol^U. lU j 


OZo 


n /I a I'D r\«\ 
U.4ol^U.U0 ) 


DDo.oo 


OO 1 




U.O 1 ^U.U4 ) 


1 OQ/1 
iZo4 


u.yy (^u. lu ; 


iZ /O 


U.OOl^U.U / ; 


( D. ( y 


oD.D 




U.40^U.U0 J 


iOio 


U. M I^U.Uo ) 


QQA 


n OQf n r\A \ 

U.Zol^U.U4 ) 


DoO. it) 


40. U 




n /1 1 /'n f\A ^ 

U.4H^U.U4 J 


iocSO 


U. i l^U.U ( ^ 




n 1 cc/'A no^ 
U.lOl^U.UZ ^ 


Dyo.D / 


o / .z 




U.4ol^U.U4 ) 




n Qi /"n 1 
u.yii^u.iu J 




U. iol^U.Uo ) 


oyy.D / 


Oo. z 




Q/i/n r\A\ 

U.041^U.U4 ) 


Z U O 


U. / Dl^U.Uo ) 




A 1 /I /A 


f UU. ( o 


oy.z 




U.oD^U.U4 J 


looU 


n 77/0 n2\ 
U. M I^U.Uo ) 


DZU 


n 1 fi/'n n9\ 

U. IDI^U.UZ j 


702.73 


61.0 




0.41(0.04) 


1978 


0.67(0.07) 




0.11(0.02) 


712.73 


70.5 




0.41(0.04) 


2357 


1.02(0.11) 




0.11(0.02) 


722.71 


79.9 




0.34(0.04) 


2413 


1.03(0.11) 


490 


0.12(0.03) 


724.66 


81.7 




0.34(0.04) 


2137 


0.99(0.10) 


1067 


0.19(0.03) 


725.65 


82.6 




0.34(0.03) 


2260 


1.02(0.11) 


568 


0.11(0.02) 


726.66 


83.6 




0.35(0.04) 


2322 


1.10(0.11) 




0.08(0.02) 


727.67 


84.5 


160 


0.32(0.03) 


2364 


1.02(0.10) 




0.10(0.02) 


728.67 


85.5 




0.20(0.02) 


1802 


1.24(0.14) 


1127 


0.19(0.03) 


729.67 


86.4 




0.37(0.04) 


2687 


1.15(0.12) 


680 


0.14(0.02) 


737.70 


94.0 




0.29(0.03) 


1941 


0.85(0.09) 


459 


0.09(0.01) 


755.62 


110.9 




0.32(0.03) 


2236 


1.14(0.12) 


1031 


0.16(0.02) 


799.52 


152.2 


525 


0.51(0.05) 


3809 


2.18(0.22) 


1669 


0.23(0.03) 



^FWHM is not presented when the spectral resolution is bigger than the measured value. 
^la uncertainties are given in parentheses. 
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Fig. 1. — r'-band image (3.5' x 3.5') of the field around SN 2005gj obtained with the Swope- 
Im telescope at LCO. North is up and east is to the left. Sixteen comparison stars in common 
between SDSS and CSP used to derive differential photometry of the SN are labeled as in 
Table [H The SN is close to the center of the field. 



T 




— Ca II Fe III 



I I I I I I I I I I I I u 

4000 6000 8000 

Rest wave length (A) 

Fig. 2. — Spectra of SN 2005gj obtained from Oct. 1 (~7 days after explosion) to Nov. 28 
(~61 days after explosion) of 2005. The sequence show the dramatic spectral evolution of 
the SN from a very blue continuum with strong Hydrogen-Balmer lines in emission in the 
early phases, resembling the spectrum of a Type Iln SN, to a Type la supernova-dominated 
continuum with broad absorption and emission features (P-cygni profiles) of blended Fe II 
and Fe III profiles. The spectra are shown in logarithmic flux scale and a constant shift 
has been apphed for clarity. The wavelength is in the rest-frame corrected using z — 0.0616 
for the host galaxy. We show the UT date when the spectra were obtained and the epoch 
(rest-frame days after explosion) in parenthesis. 
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T 




Ca II/O I 



I I I I I I I I I I I I I I I 

4000 6000 8000 10^ 

Rest wave length (A) 

Fig. 3. — Late time spectra of SN 2005gj obtained between Dec. 8, 2005 (~71 days after 
explosion) and Mar. 6, 2006 (~152 days after explosion). The labels, axis and symbols are 
the same as in Figure [2j The earth symbol shows the position of a telluric feature present in 
some of the spectra. 



-51 - 



=F I r 



n I I r 



n I I r 



n I I r 



n I I r 



10 



15 



CO 



GO 

o 
o 

+ 

Of) 
CO 

a 20 





4k 
















Am 








A^ 



OOSECPCP® O Q 



X A 



o 

o 

o o 
o 

°o 
o 

°o 
o 
o 



o □ 



[3) BSSSCQCtlQ, . 

A ^ 

ED ocaucoo- AA a 



o 
o 



O °Q 



^ AA 



°'^A 

A 

o SDSS ugriz * * 

A CSP u'g'r'i'YJHK^ 
□ MDM griz 



K-6 



H-5 
J-4 

Y-3 
z-2 
i-1 
r 

g+1 
u + 2 



25 



U I L 



J I I L 



J I I L 



J I I L 



J I I L 



650 700 750 800 

JD - 2,453.000 (days) 



850 



Fig. 4. — Observed light curves of SN 2005gj from SDSS (open circles), MDM [open squares) 
and CSP/Swope {filled triangles). The error bars are smaller than the symbols. For clarity, 
the light curves have been shifted by an arbitrary constant. 
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Fig. 5. — Difference between the synthetic g — r color calculated from the spectra and the 
observed color from the photometry. We do not include the latest spectra obtained on Jan 
24. and Mar. 6 because there is no contemporaneous photometric data. 
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Fig. 6. — Time evolution of the colors of SN 2005gj {filled circles). For comparison we also 
show the color evolution of the overluminous Type la SN 1991T (solid line), the Type Iln 
SN 1999el (dashed line), and two previous cases of Type la strongly interacting with its 
circumstellar medium, SN 2002ic (stars) and SN 1997cy (triangles). 
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Fig. 7. — Absolute ugri light curves of SN 2005gj {filled circles). For comparison we also 
show the absolute hght curves of the overluminous Type la SN 1991T {solid line), SN 2002ic 
{stars) and SN 1997cy {triangles). The error bar in the lower right pannel represents the 
typical error in the absolute magnitudes dominated by a 10% uncertainty in the Hubble 
constant. 
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Fig. 8. — Absolute light curves of SN 2005gj in the NIR: J {top panel), H (middle panel) 
and Kg [bottom panel). For comparison we also show the absolute light curves of a normal 
Type la [solid line) and the Type lln SN 1999el [open triangles). The error bar in the lower 
right of each panel represents the typical error in the absolute magnitudes dominated by a 
10% uncertainty in the Hubble constant. 
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Fig. 9. — Examples of black-body fits {solid line) to the SED of SN 2005gj obtained by trans- 
forming the rest-frame ugriz magnitudes to monochromatic fluxes at the effective wavelength 
of the filters {filled circles). These examples show the quality (i.e., goodness-of-fit) range of 
the black body fits at different epochs: =0.3 {left panel), 2.4 {middle), 4.4 {right). The 
units of flux density in the y-axis are mJy = 10^^^ ergs~^ cm~^ Hz~^. 
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Fig. 10. — Top panel: Quasi-bolometric {open circles) and bolometric light curves of 
SN 2005gj {filled circles). The bolometric light luminosities were obtained after applying 
bolometric corrections calculated from black-body fits to the optical SED obtained from the 
ugriz photometry. The dashed line shows the best-fit linear decay of 0.014 magday"^. The 
middle and bottom panels show the evolution of the black-body temperature and radius, 
respectively. 
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Fig. 11. — Comparison of spectra of SN 2005gj at 26, 46, 62 and 83 days after explosion 
with comparable epoch spectra of SN 2002ic (Hamuy et al. 2003) and SN 1997cy (from 
SUSPECT database). The spectra are plotted on a logarithmic flux scale and shifted by an 
arbitrary constant. The wavelength was shifted to the restframe using z = 0.0616 of the 
host galaxy (Aldering et al. 2006). 
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Fig. 12.— Results of SNID. We show the spectra of SN 2005gj at four epochs (26, 37, 61, and 
85 days after explosion) and their best three cross-correlation library spectra. The spectra 
are plotted on a logarithmic flux scale and shifted by an arbitrary constant. 
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Fig. 13. — Fits to the spectra of SN 2005gj. We model the spectra as the sum of two 
components: (1) SN 1991T spectrum at the same epoch after explosion as SN 2005gj scaled 
by an arbitrary constant {blue line); (2) fourth order polynomial {green line). The results of 
the fits are in red and the spectra of SN 2005gj, corrected by Galactic extinction in the line 
of sight are in black. The epochs of the spectra are shown in the upper right of each panel. 
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Fig. 14. — Results from the Gaussian fits to Ha and H/5 emission features as a function 
of time. Top left panel: FWHM of the Ha-broad and H/3. Top right panel: Luminosity of 
the narrow and broad Gaussian components of Ha. Bottom left panel: Luminosity of H/3. 
Bottom right panel: Balmer decrement, ratio of total fluxes in Ha and H/3 lines. 
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Fig. 15. — Line profiles of HjS (top) and Ha (bottom) in the two highest resolution spectra 
of SN 2005gj obtained with WHT+ISIS on day 44 (left) and Magellan+LDSS-3 on day 
115 (right). The features show clear P-cygni profiles with weak absorption minima at ~ 
—200 kms~^, demonstrating the presence of a slowly moving outfiow. 
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Fig. 16. — Identification of lines in tlie spectrum of SN 2005gj obtained at 17 days after 
explosion (2 days before g maximum). The red line shows the best fit synthetic spectrum 
generated with the SYNOW code. The lines of SN 2005gj are typical of SN 1991T around 
maximum light {blue line), and very similar to the spectrum of SN 2002ic around maximum 
{green line). All the spectra have been locally normalized. We have subtracted a constant 
value to the spectra of SN 2002ic and SN 1991T for clarity. 



